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ABSTRACT 


The role of strong atmospheric forcing events in determining the 
evolution of the upper ocean during the fall and early winter cooling 
season was investigated. The historical series of surface and near- 
Surface marine observations at three mid-latitude ocean weather ships 
LPAPA (OWS P), NOVEMBER (OWS N), and VICTOR (OWS V)] support the hypo- 
thesis that the integrated effects of these events dominate This 
evolution. For example, periods when the mechanical forcing was 
greater than the long-term mean accounted for approximately 35% of 
the time in the record examined at the three stations. However 85%/ 
68%/57% of the sea-surface temperature change at OWS N/OWS P/OWS V 
occurred ‘during These periods. 

Forty-nine data sets were examined and modeled during periods of 
intense fall and. winter forcing. The significant thermal structure 
modifications observed during these strong events were simulaved 
successfully using three modifications of the Kraus and Turner (1967) 
One-dimensional model. Evidence is presented which demonstrates That 
the amount of mechanical ly-generated turbulent kinetic energy avail- 
able for entrainment decreases as the mixed-layer depth increases. 
Furthermore, in agreement with Gill and Turner (1976), these case 
studies suggest that only a small percentage of the convectTively- 
generated turbulent kinetic energy is available for increasing the 


potential energy of the ocean by entrainment. 
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|. INTRODUCTION 


peeeerURPOSE OF THE STUDY 

The fundamental objective of this study was to investigate the role 
of strong atmospheric forcing events In the modification of the upper 
ocean thermal! structure during the fall and early winter cooling seasons. 
Simpson (1969) has demonstrated that the Significant air-sea exchanges 
(heat, moisture, and momentum) in mid-latitudes are concentrated almost 
entirely into synoptic-scale forcing events. For example, during a three- 
month period of strong winter storms at ocean weather station CHARLIE 
MmeeeoN, 355.9W) in are Atlantic, 84% of the evaporation took place in 
only 30% of the time intervals. The mechanical energy and sensible heat 
exchange were similarly concentrated. Additionally, these large forcing 
events were identified with the travelling extratropical cyclone families 
that dominate the mid-latitude weather maps over ocean regions. 

Simpson's analysis was directed, however, toward understanding the 
role of the large heat fluxes in modifying the cyclones themselves. The 
response of the upper ocean to these large heat and energy fluxes was 
not considered. In fact information, found in the literature, regarding 
The formation and destruction of transient tnermoclines by these large 
atmospheric forcing events is rather qualitative in nature. A detailed 
investigation to determine the significance of these large storms to 
the total evolution of the upper ocean thermal structure is therefore 


needed. It was the purpose of this study to partially fulfill this need. 
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B. FUNDAMENTAL HYPOTHESIS 

Denman and Miyake (19/73) investigated the response of the upper 
ocean during a |2-day period (13-24 June 1970) at ocean weather station 
PAPA (50N, 145W) characterized by the passage of several summer storms. 
They demonstrated that the mixed layer response correlated with these 
storms was significantly larger than the changes in the thermal struc- 
ture observed before and after the storms. The success with which 
Denman's (1973) numerical version of the Kraus and Turner (1967) mixed- 
layer model predicted the upper ocean response during this period demon- 
Strated that the response was largely one-dimensional. That is, the 
heat budget of the mixed layer was closely determined by the vertical 
heat fluxes only. 

The heat and mechanical energy fluxes reported by Simpson (1|969) 
were significantly larger than the fluxes observed by Denman and Miyake 
(1973). Therefore, the fundamental hypothesis underlying this research 
is that significant upper ocean therma! structure modifications take 
place during strong atmospheric forcing events in the fall and early 
winter. Furthermore, these responses are largely one-dimensional and 
are principally the result of mechanical mixing and convective adjust- 
ment of the upper layers. We may therefore apply one-dimensional model- 
ing techniques to examine the relative importance of the vertical mixing 
and convective processes during these strong fal! and early winter 


events. 


Meee ECIFIC OBJECTIVES OF THE STUDY 
The specific objectives of this study were to examine the response 
of the upper ocean during a large number of strong atmospheric forcing 


events during the fall and early winter cooling season and attempt to: 
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(|) determine the nature of the structural! modifications that occur 
during these strong forcing events; 


(2) determine the extent that a modified version of the Kraus and 
Turner (1967) one-dimensional mixed-layer model is capable of 
simulating upper ocean response during these events; 

(3) establish the principal physical mechanisms that cause these 
modifications, and quantify the relative importance of mechani- 
cal mixing and convection; 


(4) establish the significant characteristics of the atmospheric 
forcing during these events; 


(5) determine what percentage of the total seasonal mixed layer 
response may be explained during these strong events. 


These objectives represent a substantial departure from previous research 
designed to understand the modifications that take place in the upper 
Ocean thermal structure. If is the first study designed to demonstrate 
that strong atmospheric forcing events dominate the fall and early winter 
erosion of the thermocline. Additionally, it is The first research that 
attempts to simulate the response of the upper ocean during these strong 
fall and winter events. 

The formulation of a new one-dimensional mixed layer model is not 
an objective of this thesis. To a large extent the modern one-dimensiona! 
theories, developed during the past decade, have had only a few isolated 
tests with real data. This research will therefore be directed toward 
validating a number of modern mixed layer models with a large number of 
data sets gathered at the Three North Pacific ocean weather stations. 
It will be demonstrated that the Kraus-Turner (1967) one-dimensiona! 
theory may be modified to adequately simulate these strong fal! and winter 
forcing events. Additionally, the modified Kraus-Turner mode! will be 
used To isolate and examine the relative importance of mechanical mixing 


and convection at the three ocean weather stations. 
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Dee OVERVIEW OF THE THESIS 

In the next chapter the fundamental principles governing the evolu- 
tion of the upper ocean will be reviewed and the important assumptions 
used in this thesis examined in detail. This chapter will further serve 
as a review of recent mixed-layer modeling theories. The ocean areas 
studied in this research included the regions occupied by ocean weather 
Ships PAPA (50N, |145W), NOVEMBER (30N, |140W), and VICTOR (34N, |164E) in 
the North Pacific Ocean. In Chapter ||| the general characteristics of 
the atmospheric forcing and upper ocean thermal structure found at these 
three locations will be described. The analysis techniques, employed in 
This thesis, will be examined in Chapters IV and V. In Chapter IV 
several mixed-layer models will be employed to examine a large number 
of data sets in an attempt to accomplish the first three objectives. 
In Chapter V information relating to the final two objectives will be 
extracted from the historical surface and near-surface marine observations 
at the three ocean weather stations using a new analysis technique. 
Finally, in Chapter VI, the significant findings of this research wil | 
be presented and discussed relative to their importance in understanding 


and predicting upper ocean thermal! structure modifications. 
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ll. REVIEW OF MIXED-LAYER MODELING THEORIES 


A. GOVERNING EQUATIONS 

The changes that occur in the upper ocean are governed by the conser- 
vation laws of mass and momentum, by the equation of state for sea water, 
and the laws of thermodynamics. In this section the equations represent- 
Ing these physical laws will be examined, simplified, and transformed 
into parameterized expressions from which numerical solution may be ob- 
tained. Since many of the important objectives are based on the results 
of mixed layer models, a careful review of the mixed layer theories wil | 
be presented. 

The conservation of momentum and mass are represented by the Navier- 
Stokes equation of motion, invoking the Boussinesg approximation, and 


the condition of incompressibility. 





Bu, 3u ee au. , 
woe tule te... Nu, = —— 6... ty - — qo (2-]|) 
ot J + ikeaaeck oF a ij dy Z P, 1 
Ou, 
x = 0 (2-2) 
where ees 152,25 
and u; = (Cu) ,U5,U) = (u,v,w) 
= (X 1 X59X-) = Oye) 


Since the Boussinesq approximation assumes hydrostatic equilibrium for 
the reference state of the ocean, the pressure (p), and the density (p) 
represent departures from this state. The reference density is approxi- 


mated, with sufficient accuracy, by a specified constant value OQ, 


ed 





le 026 om/cm”). Additionally, ae 


the earth's rotation vector, St is the Kronecker delta, v is the 


is the permutation Tensor, 7 is 


kinematic viscosity, and g is local gravity. Cartesian coordinates 
will be used throughout the development and the z-axis will be taken as 
positive up from the sea surface. 

In the upper ocean the density is primarily a function of tempera- 
meet) ), and salinity (s), such that a simplified equation of state 


may be assumed: 


o(x,,t) = p{l-aLl T(x, ,t)- Vg 8Ls(x.,t) meal, (2-3) 
where 
By SED = ieee: 
a= eee and B ee 


The coefficients, a@ and 8, are Taken as constants throughout this thesis 
oem = 2.5 x lon *(°e) andes a=- 7 ox lone coe 
It will be assumed that the frictional generation of heat and mole- 
cular heat transfer processes are negligible compared with typical values 
of radiant solar energy and the turbulent heat fluxes exchanged between 


The ocean and atmosphere. Therefore a simplified form of the first law 


of thermodynamics may be expressed as: 


dT iets ll din) 


=e tu (2-4) 


j ox. OC. oz 
J j oo 9 
According to Jerlov (1968), the total downward irradiance in the spectral 


range of 300-2500 nm decreases to 50% of its surface value in the first 


meter of the ocean, virtually irrespective of the water type considered. 


ie 


Therefore, in this study, the absorption of short wave radiation below 


the initial meter will be approximated by: 


R(z) +R, ef (2-5) 
where y represents a total extinction coefficient (taken as 0.4 roe 
nn represents the surface absorption, Ae =e oe » where We is the 
total solar flux at the surface. 

The final expression necessary to complete the system is the equa- 
tion for conservation of salt. Once again neglecting molecular diffu- 


Sion, if may be simply represented as: 


= tu. =0 (2-6) 


B. THE ONE-DIMENSIONAL HYPOTHESIS 

The basic assumption of the one-dimensional hypothesis is That the 
ocean is horizontally homogeneous in all its properties Cu. , Si) dione 
This assumption restricts the domain of applicability to time and space 
scales over which the vertical fluxes of mass and momentum dominate The 
horizontal fluxes. Because of the sparsity of open ocean measurements 
This assumption is difficult to verify a priori. However, the one-dimen- 
sional hypothesis is desirable at this point because three-dimensional 
models are not only far more complicated and expensive, but the frequency 
Of observations (especially oceanic) prohibits proper initialization, 


calibration, and validation. 
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ee TURBULENT FORM OF THE BASIC EQUATIONS 

The turbulent components are introduced into the basic equations 
through the Reynolds decomposition technique, whereby the variables 
are expressed as a time averaged mean and a fluctuation about the mean 


(i.e. T=T+7T'). The prime denotes the fluctuation and the overbar 


represents The mean; 


= 


_ _ 
ee 7 dt , (2-7) 
At 
2 


for example. The integral time scale, At, should be short compared 
with the time in which the mean field properties are changing, but long 
relative to the time scale of the fluctuations. 

Applying this technique to the basic set of equations (neglecting 


the mean vertical motion W ) results in the fol lowing: 











= = oe 
<< + ife = - ee (Pas 
oe | 8 a : 
| pC. az (R(z) or Wea) (2-9) 
OP 
' t 
oe Ne a (2-10) 


Equation (2-8) is the equation of mean motion expressed in complex 
notation (c = u + iv) with the geostrophic component removed. Equa- 
tions (2-9) and (2-10) are the conservation relations for mean tempera- 
ture (T) and salinity (s), while expressions wic', w'l', and wis! 
represent the turbulent fluxes of momentum, temperature, and salt re- 


spectively. The parameterization of these fluxes will be accomplished, 
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eran, bY examining the nurbulent kinetic energy budget for The upper 
ocean, which is derived from the Navier-Stokes equation. 

lf (2-1) is multiplied by u, and subjected to Reynolds decomposi- 
tion, the resulting expression wil! represent the conservation of total 
kinetic energy. The conservation of mean kinetic energy is formulated 
by the decomposition of (2-1) multiplied by the time averaged mean 
velocity U, . The turbulent kinetic energy equation is obtained by 


subtracting the mean from the total equation and the resulting expression 








er 
2 2 dU 
9g Opie a na see Z 
5 (>) eet G ts 5 7) oe exe ; € CZ =|") 
O O 
where a = utu.l, 
Ue = UN Ue 
a 
f 
and c= u(x )? 
Oo 
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The first term on the left is the local time rate of change of Turbu- 
lent kinetic energy while the second, the divergence term, specifies the 
vertical redistribution of the turbulent kinetic energy wi(pt/p + a2 
by the turbulence. This energy is generated primarily by breaking waves 
at the Piicface and by the Reynolds stresses acting on the mean flow. The 
right hand terms represent the three ways in which turbulent kinetic 
energy is either gained or lost. The first of these, which is usually 
positive, is the rate at which mean kinetic energy is converted to turbu- 
lent kinetic energy by the working of the Reynolds stresses “2, u_twt 


3U me 


against the mean velocity gradient = - The second, the covariance 


oe 





between the fluctuations tn nee and the vertical velocity, may be 
positive or negative. If the basic density distribution is statically 
unstable, then fluid aiencets moving upwards tend to be less dense than 
those descending, and a release of potential energy takes place by free 
convection (-g wrpt/p. > 0). %If, on the other hand, the density distri- 
bution is stable, the reverse is true and -g TS oe vee inis covariance 
fem represents the rate at which Turbulent kinetic energy is expended 

5 mixing the less dense fluid downward, thus increasing the potential 
energy. The last term (€) is the rate at which turbulent kinetic energy 
is dissipated by viscosity and always represents a loss. To solve the 
system (2-8) - (2-11) requires specification of the vertical structure 

of the properties and the boundary conditions imposed on the domain. 
Moreover, this involves parameterization of turbulent processes and addi- 


Tional simplifications will be necessary to make the problem tractable. 


freee BULK MODEL HYPOTHESIS 

In this thesis the one-dimensional dynamic processes that affect the 
evolution of the upper ocean will be studied in terms of the energetics 
associated with the turbulent kinetic budget. The parameterization of 
These energetics is accomplished by idealizing the upper ocean structure 
with the assumptions of the bulk model hypothesis (depicted in Fig. 2-1). 
The quantities Co, Ts, Se represent the vertically averaged values of 


velocity, temperature, and salinity, defined for example as 


0 


gas x 
ce = eee 7 Cedz (2-12) 
-h-6 


With this concept, the density structure (To, S.) directly below the 


wind blown ocean surface is assumed vertically homogeneous to a depth of 
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Figure 2-I|. Schematic of the oceanic mixed layer. 
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z=-h. Below this level there is a discontinuity in the density fol- 
lowed by a stable density profile. The mean velocity structure is 
modeled as being vertically uniform and fully turbulent in the mixed 
layer and negligibly smal! and nonturbulent below. Deviations from this 
vertical structure occur at the top and bottom of the layer where shear 
zones are formed. These shear zones are formed at the top due to the 
action of the wind on the water, and at the bottom due to the slab-like 
motion of the mixed layer over the quiescent water in the pycnocline. 

The surface shear zone is known as the production zone, where turbulent 
kinetic energy is generated by the working of the Reynolds stresses 
against the mean velocity gradient. The bottom shear zone is an entrain- 
ment zone of thickness 6. There is a vertical flux of mass and momentum 
at the top (z = -h) while none leaks out the bottom (z = -h-0). It is 
further assumed that the Reynolds stresses are also capable of generat- 
ing turbulence in this zone. Thus the bulk model concept assumes that 
Variations from the structure of the upper ocean depicted in Fig. 2-1 

may be neglected and The vertical mean quantities may be predicted, with 
sufficient accuracy, by specification of the turbulent transfer processes 
at the boundaries. : 

To maintain these homogeneous profiles throughout the mixed layer a 
continual vertical flux of turbulent energy is necessary throughout the 
layer. When there is a convergence of turbulent energy at z = -h, the 
entrainment zone is destabilized and the excess turbulent kinetic energy 
is expended by entraining fluid from below as the layer deepens. With a 
downward buoyancy flux at the surface (excess heating or precipitation), 
it is possible that an insufficient flux of turbulent kinetic energy may 
be available to mix the fluid homogeneously to the existing mixed layer 


depth. In this case a new mixed layer depth is established at the level 
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where the downward vertical turbulent flux vanishes. Since this level 

is higher than the previous layer depth, this formation is called layer 
retreat. The turbulent motions below this level are assumed to be shut 
off from the energy source, and become nonturbulent by viscous forces 

on a dissipation time scale. IT is therefore envisioned that the mixed 
layer is being constantly re-established from the surface, and that 
Shallowing mixed layers are not the result of the interface moving upward, 
but are ene consequence of net surface heating and an insufficient down- 


ward flux of turbulent kinetic energy. 


E. BOUNDARY CONDITIONS 
The boundary conditions will be specified as a function of time and 
the overbar on the mean quantities will be dropped. The surface stress 


is denoted as cs and 


ae iO 
- w'c'(o) = (2=15) 
Po 
Onecare 
while ics Wop aaco See (2-|4) 
pac ae 
Op Op 


The term Q. is The sum of the turbulent fluxes of latent (Q,) and 
sensible (Q,? heat and The effective long wave back radiation (O.). 
The short wave solar energy (Q.) is always taken as a positive value. 


The surface flux of salt is 
- w's!(o) = s(E-P) (7=15) 


where P is the rate of precipitation, and &€ is the evaporation rate. 


Zo 


The boundary conditions at the bottom of the mixed layer are derived 


Baeeaiegrating (2-3) - (2-10) over the entrainment zone (-h-d < z < =h), 


Thus 
—— dh 
tanto = =< — 
Wie (iy y+ Ac 
where Nes = CC) melo.) 
Since only one-dimensional effects are considered in this model, oh 


is the time rate of change of the mixed layer depth (h) due to turbulent 
processes. As this expression implies an upward momentum flux for 
dn 


i. Peecan impossibility in this system), if is rewritten as 


wic'(-h) = - A = Ac (2-16) 


where A is the Heaviside unit step function defined as 


Oars ink on < 0 
A = (2-17) 
Ae \eb = 7 


The fluxes of heat and salt at tne base of the mixed layer are 


Similarly derived and are 





R es 
MEO art | (ony Ss oe ey OM AT (2-18) 
p o_C ot 
O-D O-D 
i dh 
f | Gar eae - . een = 
wret(-h) = -A Spd (2-19) 
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F, THE INTEGRATED MASS, MOMENTUM, AND ENERGY BUDGETS 
Equations (2-8) ~ (2-10) may now be integrated over the mixed layer 


to form the mass and momentum budgets of the region. The results are: 


as eoleee 2 ah 
ace Bi ea aoe, C2=720)) 
-0 -R @Y" 
Sis 1 [Ss Nao? ae (2=21) 
ot h pC ot 
OP 
dS 
an re ; 
Ae = i s« P) -A y+ ss (2-22) 


These equations are subject to simple interpretation. Equation (2-20) 
shows that the mean motion of the mixed layer is modified in time by rotation 
and the vertical fluxes of momentum at the surface and base of the slab. 
The surface flux may add or subtract momentum depending on the relative 
directions of T. and C. . During the period when the mixed layer is 
deepening the flux of momentum at the base is always a sink for energy 
(Ac = C. io) and represents the cnergy necessary to impulsively “aecelc- 
rate the entrained fluid to the velocity of the mixed layer. For deepen- 
ing mixed layers the density, specified by (2-21) and (2-22), wil! 
normally increase due to the turbulent flux of density at the surface 
and base of the mixed layer, while it normally decreases for shal lowing 
layers. For stable temperature and salinity structures, deepening layers 
normally become cooler and more saline while the reverse is true for the 
retreating case, 


To close this system the time rate of change of the mixed layer 


depth (dh/dt) must be specified. This may be accomplished through the 
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integration of the turbulent kinetic energy equation over the mixed 
layer. It is at this point where most modern one-dimensional bulk 
models differ and the procedure will be given careful consideration. 


One may start by defining the terms of (2-I1) as follows: 


. 9 eae aU, 
aa eat a 2 a dt = 
Gy, = -2, I& Cw ee )] + Uw? a} az 25) 
-h-6 
O 
= P, f E dz (2-24) 
-h-6 
O 
eS = P, w'b? dz (2-25) 
-h=6 
0 ay 
0. oe) dz (2-26) 
=h-6§ 


itemminsi term (G,) represents the total contribution to the Turbulent 
kinetic energy budget by mechanical production processes. The second 
expression (D,) is the total dissipation in the layer by viscous forces. 


The rate at which potential and turbulent kinetic energy are being ex- 








eyenged is represented by P, , where w'b' is the turbulent flux of 
buoyancy. 
wio! aa a 
w'b! = =-g = gtd Wit’ = 8B aw's*) (2-27) 
O 


Finally, Sy, is the rate at which the turbulent kinetic energy budget 
for the layer is changing, and is commonly known as the storage term. 


The integrated form of (2-11) may be simply expressed as 


P, + G, - Dy - Sy, = 0 28) 


28 





Kraus and Turner (1967) were the first to formulate a mixed-layer 
model based upon (2-28), except they neglected the storage term and 
assumed the density a function of temperature only. In deriving an ex- 
Maecoicon for FP, , the assumption that temperature and salinity remain 
vertically homogeneous throughout the mixed layer, and The relation- 
Ships represented by (2-9) and (2-10), require that the vertical fluxes 
of heat and salt be a linear function of depth. With this restriction 
The trapezoidal rule may be applied to these fluxes and an expression 


fOpeer, derived as 





O O 
P, = pga f w'T! dz - 0.08 ff tat 
-h -h 
; pygh Cae somes cme 
Z OG ot 
Op 
Boh 
a O hAb , oh 
+ = = [s(e- P)+A— ; As] | = =a = A xe , (2-29) 
where 
- Q. - R_F(yh) 
_ S ‘ae o B 
BO = oo pa i s(e-P)| , (2-30) 
O Pp 
Ab = PgaAT = p98 NSS (2-31) 
and 
eh) = = igs, Cae aa aD) 
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The process modeled by a A on 


kinetic energy expended to deepen the mixed layer and increase the poten- 


represents the amount of turbulent 


Tial energy of the column of water (raise the center of gravity) by en- 
Bh 
Training denser fluid from below. —— jis the turbulent kinetic energy 


Z 
(released potential energy) generated when BS < 0 and free convection 
occurs. When BS > 0, it represents the energy expended during forced 
convection to mix the buoyant surface water downward and increase the po- 
tential energy. 

Tne function F(yh) reflects the erence of solar energy to some 
depth with the property that F(yh) > 0 as yh+o (complete absorption 
In the layer), and Fl(yh) > | as yh > O (complete penetration through 
miemiayer). Substituting (2-29) into (2-28) results in an expression iow 
The Time rate of change of the mixed-layer depth in terms of G, , Dy, , 

S 


feeena ane surface tluxes of heat and salT; 


A we = + =(E-P) + 
ot (AT- Eas) 85s Oo pgah 


” : ae 3 | a. 
When the right hand side is positive a deepening mixed layer is predicted. 
When the terms are negative the Heaviside function is equal to zero and 
a diagnostic equation is formed to obtain the depth to which the layer 
Peireats. The methods by which 5S, , G, , and D, are parameterized wil! 
now be described. 

The storage term is parameterized in a manner described by Kim (1976). 
The vertical mean turbulent kinetic energy of the layer, postulated as 


being quasi-invariant in time, is defined as 
C 2s i. qé dz (2~34) 
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and the storage term parameterized as 


ss 2 on 
oe aemee C + (2-35) 
Kim further assumes that 
C~ =m, wy. ese) 
m | 
where m is a nondimensional constant, and 
ae = b3) wea 
Po 


is the friction velocity of the water. 
The parameterization of the total mechanical production of turbulent 
kinetic energy may be accomplished by integration of (2-23) over The 


mixed layer, with the integration intervals depicted in fig. 2-I. 


=i 


3U 
Gy = -9 Ae +¥) (e * FI + Lr reas f ao | 2-5) 


~h=6 


The first two terms represent the mechanical] production at the surface 

due to atmospheric perturbations and breaking waves (term |), and The 

work performed against the mean velocity gradient by the Reynolds 

stresses (term 2). The third term is the shear production in the entrain- 
ment zone as the layer deepens. Niiler (1975) carefully performed these 


integrations to show that the contribution of the first two terms is 


Sal 


> 
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Seopertional to the cube of the friction velocity (w,), while the third 


; dh IC. 
is approximately A 35 son ae Therefore 
Col 
_ iS oh | S 
ee = M5PGW x a5 mz0 A ay oe (2-39) 


and the nondimensional constants m5 and mz are both of order unity. 
The total integrated turbulent kinetic budget is obtained by substi- 


mma (2-29), (2-35), and (2-39) into (2-28). 
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poe ge a A~=- oo Cr De (2-40) 


G. CLASSES OF MIXED LAYER MODELS 
|. The Prototype Turbulent Bulk Model 
Most modern mixed layer models can - classified as to which 
terms of (2-40) are used to predict changes in the mixed layer proper- 
Ties. Kraus and Turner (1967), hereafter KT, formulated the prototype 
Turublent bulk model, assuming a balance between terms ||, II11, and V 


in (2-40) and neglecting the rest as follows: 


ea 
hAb oh _ iS O 
we og oS ze ae 
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When the sum of the terms on the right are positive, the layer deepens 


according To 


3 
an _ Pot 7 To” (2-42) 
at hb 


However, during periods of weak winds and surface heating, the right 
hand side becomes negative and 
20 Wy 
ho == (2-43) 
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is a prognostic equation which calculates the retreating mixed layer 
depth proportional to the Monin-OQbukhov length scale (L). Incidentally, 
Kitaigorodski (1960), using a steady state model and dimensional! analy- 
Sis, reasoned that this should be ine proper length scale. 

Kraus and Turner demonstrated that this model was capable of 
Simulating the annual evolution of the mixed layer with a saw tooth 
heating function and constant wind stress. Denman and Miyake (1973) 
applied a numerical version of KT to a |2=<day period at ocean station 
PAPA with favorable results. Using observed forcing they were able to 
simulate both the daily and weekly changes in mixed layer temperature 
and depth during a period of moderate synoptic scale winds, apparently 
without any significant effects from vertical or horizontal advection. 
On the other hand, Dorman (1974) applied the same model to cases of 
spring heating and fall cooling at ocean station NOVEMBER with | imited 
Success, and attributed the model's poor performance to horizontal! 
advection. 

However, these few simulations are not conclusive evidence for 


The validity of the one-dimensional bulk model (in particular KT) and 


She, 
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a primary objective of this work will be to attempt a large number of 
experiments to gather such evidence. In addition to KT, the model of 
Kim (1976) and the model of Elsberry, Fraim, and Trapnel! (1976), here- 
after KIM and EFT, will be evaluated. Both of these models assume the 
dominance of surface production and are therefore of the KT Type. 

Kraus and Turner originally determined that KT predicted exces- 
Sive mixed layer depths when mechanical production and free convection 
were considered together. They concluded that either their parameteri- 
zation of mechanical production was in error or the effects of dissipa- 
tion could not be neglected, or both. Turner (1969) found that the 
fraction of the total energy, imparted during an impulsive wind event, 
That is used to increase the potential energy by entrainment was larger 
Than OW, , as originally reasoned by Kraus and Turner from dimensional 
analysis. He suggested that much of the kinetic energy goes into drift 
currents and is eventually used to deepen the layer; however he did 
not specify the mechanism for this deepening. 

Beene Prototype Inertial Bulk Model 

In an attempt to explain the rapid deepening of the mixed layer 
in response to strong impulsive forcing, as reported by Turner (1969), 
Pollard, Rnines, and Thompson (1973), hereafter PRT, formulated a model 
quite different from KT. Assuming the density was a function of tempera- 
Ture only, they neglected the turbulent kinetic energy budget and con- 
sidered only the time rate of change of total kinetic (KE) and potential 


(PE) energy as follows: 


az t= = Tuto) (2-44) 


34 





‘7 


11 


é 


Peni 


= irene. 


7 





Further, they implied that at the onset of heavy winds the mixed layer 
would respond and move as a slab, and through the mechanism of mean flow 
instability the mixed layer would deepen. They postulated that the mean 
flow would remain unstable as long as the bulk Richardson number (Ro) 


was less than or equal to unity. 


* _ gahaT zy 
Co 
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R, 


. | (2-45) 


The PRT model predicts continual deepening as long as Tuo) is posti- 





% 
tive and R, <!l. Freecto)= ByeriS > 2 Pari VeUlcn «SO ul Kenmare 
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At time t=mT/f (half the inertial period), TUCO) becomes negative, 
the energy flow to increase h ceases, and since the water cannot unmix, 
* 
h must remain constant and R > | . The PRT model predicts a maximum 
mixed layer depth for a constant w, as 
Reais, 


h —— (2-47) 
max JEN 


where N* = gal and N is the Brunt-Vaisala frequency. 
Niiler (1975) concluded that the PRT model, in reality, assumes 

a balance between terms !I| and !V in (2-40). If ne= | , then this 

balance simply reduces to (2-45). The model is most effective in causing 

mixed-layer depth changes during periods of relatively shallow mixed 


layers and strong impulsive winds. However,a major difficulty with this 


a5 


model is that it fails to account for the gradual deepening that takes 
place during periods of weak forcing or when The initial layer depth 
is greater Than apes . It is entirely possible that the deepening of 
the mixed layer may occur as a combination of KT and PRT and that the 
geophysical situation (magnitude of the wind, layer depth, and ocean 
Stability) will dictate which process will dominate. 

In section IV the one-dimensional bulk models of KT, EFT, and 
KIM will be applied to data sets obtained during the autumn and winter 
meen in the North Pacific. The PRT model was not included in this 
study because the calculations shown in Table 2-| indicate that the bulk 
Richardson number remains greater than unity during all but a small num- 
ber of observations during these periods Che > 50 m2) In Ser tonming 
these calculations it was assumed that at time +t = 0, h = ii ,» and 


both PRT and KT type deepening occur simultaneously. Additionally sur- 


face heat fluxes were neglected, a constant temperature gradient was 


ht x 
assumed, and AT was approximated by i The values of R. (as 
Peemer ion Of wind speed) are calculated at + = 1/f . Therefore values 

* 
of R. < | indicate situations when the PRT model would have deepened 


a greater amount Than the KT model. Also, it has not been clearly 

* 
established that the mean flow positively becomes unstable at Rs = | 
and, in fact, the instability may be initiated at much lower values (see 


iiemen, 19/5, op. 97-102). 


H. NON-PENETRATIVE FREE CONVECTION 

Returning to the KT model two difficulties are yet to be resolved. 
The first of these involves the percentage of turbulent kinetic energy 
generated during free convection that is actually utilized for entrain- 


ment. Kraus and Turner followed Ball (1960) and assumed that 100 percent 
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feeee 2—-|. R as a function of wind speed (Uso) and inittal layer 


depth (ho). 
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is used. However, laboratory work by Deardorff, Willis, and Lilly 
(1969), and a more recent field experiment by Farmer (1975) indicate 
that only a small fraction of this energy (1-3%) is actually converted 
to potential energy through entrainment. Furthermore, Gill and Turner 
(1975) demonstrated that a fraction equal to 0.15 was adequate to 
achieve annual cyclic steady state in the potential energy balance. 
The reason becomes evident by considering a simple case where solar 


radiation is neglected and density is a function of temperature only. 


Then, 
O 
= = -e,g0 f w'T! dz (2-48) 
<h 
p.gah f~ Q 
or Jane O Lee ecu . 
os 5 Ar A $F ar] (2-49) 
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Sor OT = © oe * 5 gah (2-50) 


where ris the fraction of convectively-generated turbulent kinetic 


energy utilized for entrainment, (2-49) becomes 








ee ee 
oT Z oC.. * Gamat. oe (2-51) 
ae 
ieesn!| mode| assumes r= 1 and neglects DO, and 5S, . 
as = Gy = OW, (2-52) 
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Equation (2-52) demonstrates clearly that a cyclic steady state balance 
is impossible in the KT model and, except when the wind stops blowing, 
the potential energy increases continually. The consequences are exces- 
sive mixed layer depths and downward heat flux into the deep layers. 

As the entrainment process takes place, a continual downward heat 
flux into the deeper layers occurs. Before the seasonal thermocline may 
be reestablished, during the spring and summer heating cycle, this heat 
must be removed by non-penetrative free convection (cooling without deep- 
ening ). The KT model cannot accomplish this process and, when inte- 
grated over several annual cycles, will eventually erode away the thermo- 
cline. Therefore, the conclusion is drawn that non-penetrative free 
convection (r < 1) must be an integral part of any mixed layer model 
and, following Gill and Turner (19/75), a value of r=0.15 will be 
used in all the models evaluated in this thesis (including KT). The 
theory explaining how such a large percentage (85%) of the kinetic 
energy is lost is not complete. A certain amount is dissipated by vis- 
cosity, and atmospheric investigations such as Townsend (1968) 
and Stull (1975) indicate that internal waves generated by convection 
are capable of radiating energy away from the mixed layer interface 


into the gradient region. 


|. DISSIPATION ENHANCEMENT 7 

The second difficulty with KT is that excessive deepening is pre- 
dicted, even when free convection is neglected. This is a consequence 
of not properly parameterizing dissipation. A tank experiment by 
Thompson and Turner (1975) demonstrated that the entrainment rate pro- 
duced by a stirring grid is not directly related to the velocity of the 


stirrer, but to the turbulent velocity near the interface. They reasoned 
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that the turbulent kinetic energy density decays with depth and that 
less is available for mixing as the layer deepens and the entrainment 
zone gets further away from the surface production zone. 

The first model to include this process was that of Elsberry, Fraim, 
and Trapnell (1976). Total dissipation was assumed to increase 
exponentially as a function of layer depth, h, and a scale depth, 2Z, 


and 


Gy = Dy = OnWy” 54 (2-53) 


In this model G, - Dy, is nearly equal to the downward energy flux 
from the atmosphere (according to KT) If the mixed layer is shallow 
(h << Z) while it tends towards zero for very deep layers. 

A more recent model by Kim (19/76) also uses a depth dependent dis- 


Sipation parameterization and additionally includes the storage term. 


3 
Cre —. Ogee, 25 PiWe OD 25a) 
where D, 1S a constant background dissipation. The surface production 
term is parameterized according to Kato and Phillips (1969), and D. 
calculated from the dissipation data of Grant, Moilliet, and Vogel 


(1968). Kim's model has the additional interesting property of being 
able to predict steady state in the potential energy balance for neutral 


conditions (no heating or cooling). For the KIM model (2-51) becomes 
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(2-54) 
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Thus, neutral! steady state is predicted for Tt >> 





2D 
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(2-56) 


is a constant determined by initial conditions. 


Wx 
; this means 
b 


that a constant wind is only capable of deepening the layer to a depth 


described by (2-56). 


fee eR ICAL DIFFUSION 


The prototype turbulent bulk model has the undesirable tendency to 


predict abnormally large temperature gradients at the base of the mixed 


layer (see Denman and Miyake, 


\o7S,7Gi¢g. 6). Inisals a econsealenec 16; 


the zero-flux condition imposed at the base of the entrainment zone. An 


attempt wil | 


assuming that the thermocline is weakly diffusive. 


be made in. this research to overcome This difficulty by 


The temperature be- 


low the mixed layer will be specified by 
2 
a =A ow (2-57) 
a2 - 
where oy is the vertical diffusion coefficient. The KT, KIM, and EFT 
will be evaluated with this diffusion tendency to determine its effec- 
Ti veness. 
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fee OALINITY EFFECTS 

An additional assumption, made throughout this thesis, is that the 
density structure is a function of temperature only and that buoyancy 
flux (w'b!) is synonymous with heat flux (w'T') . This assumption 
is necessary (but certainly not desirable) because the data used in this 
investigation included neither salinity structure information nor ob- 
served precipitation rates. The errors introduced into the system by 
this assumption are reflected by (2-53), and (2-40). These equations 
indicate that salinity changes may be an important consideration depend- 
ing on the season, geographical location, and vertical depth scale over 
which the model is applied. An attempt will be made to determine to 
what extent this assumption affects model performance. 

In summary, the fundamental principles governing the evolution of 
the mixed layer have been reviewed in some detail, and the assumptions 
of the one-dimensional bulk model specified. The expressions that wil] 
be used to examine the physical processes responsible for observed 


changes in the mixed layer are (2-21), (2-33), and (2-40). 
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ieee Chane le mol ConOnmine  AhMOomnmen | C Gore RNG 
AND OCEANIC STRUCTURE AT OCEAN WEATHER SHIPS 
(OWS) PAPA (P), NOVEMBER (N) AND VICTOR (V) 


es 


peeeeeadA SOURCES 

One of the largest data sets for investigating the upper ocean ther- 
mal response to atmospheric forcing is the meteorological and oceano- 
Graphic observations taken at ocean weather ships. To examine the physi- 
cal mechanisms responsible for changing the thermal structure, it was 
decided to focus the analysis on the data gathered at Ocean Weather 
Ships PAPA (50N, |[45W), NOVEMBER (30N, I40W), and VICTOR (34N, 164E) in 
the North Pacific. Furthermore, since the objective of this thesis is 
To examine these processes during the fall and early winter cooling 
season, only data collected during September through December will be 
considered. 

The near-surface marine observations were provided by the National 
Weather Records Center and the three-hourly data included measurements 
of sea-surface temperature re air-tTemperature (T), dew point (TW), 
wind speed (ul), and visual estimates of total cloud cover (C). Table 
(3-|) lists the years when data were available and includes approximate- 


ly 96% of the possible three-hourly records. 


TABLE 3-1. Availability of data at the ocean weather stations. 


Station Atmospheric Observations Mechanical BT's 
F 1946, 1948-1970 1946, 1948-|970 
N 1946-195] |947-|950 
|953-|954 | 954-1970 
1956-1970 
V 1956-1970 1956-1970 
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Information regarding the evolution of the oceanic thermal structure 
was obtained through analysis of mechanical bathythermograph (MBT) records 
provided by the National Oceanographic Data Center. The years for which 
MBT data were available are also listed in Table (3-I|); however, the fre- 
quency of These observations was highly variable, and usually numbered 


less than a few hundred per season. 


B. FORMULAS FOR COMPUTING ATMOSPHERIC FORCING 

Examination of (2-21), (2-31), and (2-38) indicates that the atmos- 
pheric forcing necessary for modeling the response of the upper ocean 
includes a measure of turbulent kinetic energy flux eae the effective 
solar radiation (Q.), the turbulent fluxes of sensible (Q,) and latent 
heat (Q.)) and the net back radiation (Onde These variables may be esti- 
mated using the measured atmospheric parameters in the following bulk 


aerodynamic formulas: 


u, = von Cu, x 10%) (cm/sec) (3-1) 
Q, = 3,/6/7 Cy (On9S ET = E? u, (ly/day) (3=2) 
QO, = 2,488 C, (T= T,) uy (ly/day) (3-3) 
Q, = 1.14 x 10” /(273.16+T)9(0,39-0,05VE_) (10.607) (ly/day) (3-4) 


where u, IS the friction velocity of the atmosphere and is related to 


Wy = s? Uy (3-5) 


The non-dimensional drag coefficient (Cp) was assumed constant (1.3 x 
-3 
10 ~) throughout this study and the saturation vapor pressure of the 


marine atmosphere CE) in direct contact with the ocean was estimated 
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from the observed sea-surface temperature. The final quantity CE.) Is 
the saturation vapor pressure of the atmosphere at a height of approxi- 
mately 10 meters and was estimated from the dewpoint temperature. 
Additionally a daily estimate of the effective solar radiation was com- 
Duted from the formula developed by Seckel and Beaudry (1973). A more 
complete description of the empirical formulas is presented in Appendix 
A, along with a discussion of their WAH tac assumptions. 

C. CHARACTERISTIC ATMOSPHERIC AND OCEANIC CONDITIONS 

AT THE WEATHER STATIONS 

Before proceeding to the details of modeling the upper ocean, the 
general character of the marine atmospheric forcing and oceanic response 
at the ocean weather stations will be examined to place this study in 
perspective. 

During the fall and winter, the mean atmospheric circulation at OWS 
P and OWS V is dominated by a large barometric low pressure system 
(Aleutian Low), and the mean flow is generally westerly. The mean sur- 
face winds at OWS N are under the influence of the subtropic high and 
are generally northeasterly. However, the most outstanding feature of 
iieemid=latitude atmosphere in the North Pacific is The fluctuations in 
the atmosphere associated with the frequency and intensity of travelling 
extratropical cyclones. : 

An analysis of historical storm tracks, presented in the U.S. Navy 
Marine Climatic Atlas (1957), indicates that OWS V and OWS P lie in 
Close proximity to the preferred path of the major storm centers, while 
OWS N is located a considerable distance to the south. The influence 
that these storms have on air-sea interactions is therefore expected to 


be most pronounced at OWS V and OWS P and weakest at OWS N. 
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Figures 3-| present a comparison of the long-term mean atmospheric 
forcing and oceanic thermal response at the three ocean stations. The 
atmospheric forcing and sea-surface temperature were computed from the 
three-hourly surface and near-surface observations while the batnyther- 
mograph file was used to establish the mixed-layer depth (defined as the 
depth at which the temperature was 0.2°C less than the sea-surface tem- 
perature). Each point represents the daily average of all available 
data (see Table 3-I|) smoothed by a 7-day running mean (for display pur- 
poses only). 

These figures are presented merely to illustrate the relative magni- 
tude and variability of the atmospheric forcing and oceanic response 
that we might anticipate at the three stations. An important observa- 
tion that should be made is that the air-sea interactions in the North 
Pacific are highly dependent upon geographical location. Since we are 
interested in evaluating model performance under a variety of geophysi- 
cal situations, these three stations appear ideally suited for this 
study. It is also important to observe that there is a significant cor- 
relation between the strong, variable forcing and large oceanic response 
at OWS P and QWS V and the relatively weak, steady forcing and response 
at OWS N. In Chapter V a detailed analysis of these data will be per- 
formed, designed to understand the principal mechanisms by which these 
evolutions take place. 

The North Pacific has been divided into distinctive ocean regimes 
(Tully, 1964) based on similarities in the characteristic temperature 
and salinity structures in these regions. OWS P is located in the Paci- 
fic Subarctic region while OWS V and OWS N are located in the Pacific 
Subtropic. Figures 3-2 illustrate schematically the major features of 


the ocean structure to be found in these regions. 
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Figure 3~|. Long-term mean atmospheric forcing and oceanic 
thermal response at OWS P, OWS V, and OWS N. 


47 





/ 
¥ 3 
— ae / 
@ 
i; e 4 — 
>> Ane FAs - 
oo - 
j 
a - 
; 
i 
; f 
| 6 4 4 
j i = 


a 
. 





: ; 
- : — 
; > 
7 
mS 
J 3 é 7 
a ai 
a nt 
= OO a = Oo a a 
s eg 
Viv 
- 





TCC) SC) p(gm/cm?) 


of 6 8 10 12 32 33 34 1.024 1.026 


100 


DEPTH (M) 





300 


T (°C) Sec) p(gm/cm?) 
14 18 22 26 34 35 36 1.024 1.026 


100 


DEPTH (M) 


i) 
Oo 
Oo 





300 


Figure 3-2. Major features of the ocean thermal structure in the 
NOnmioeaeinic., (Ad Pacific Subameric Region 
(By 4racitic Subtrop ic Redon 
(after Tully, 1964). 
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The Pacific Subarctic region is distinguished by excess precipita- 
en (E-P < 0), large upward surface heat fluxes, and strong mechanical 
forcing through most of the cooling season. These quantities interact 
to create a unique temperature and salinity structure and the most dis- 
tinctive feature to be noted in Fig. 3-2A is the existence of a stable 
Salinity gradient in The upper ocean. Throughout this region the density 
structure is a function of both the temperature and salinity structure 
in the upper 100 m, with the salinity becoming increasingly more impor- 
tant at depths approaching the permanent halocline. 

At OWS N and OWS V very similar temperature and salinity structures 
are observed (Fig. 3-2B) which are very different from OWS P. The Sub- 
tropic region is characterized by large surface heat fluxes and values 
of E-P > 0 through most of the cooling season. Consequently both tem- 
perature and salinity decrease with depth in the upper ocean and the den- 
sity structure is primarily a function of the temperature structure in 
This region. 

OWS P offers a unique opportunity to examine and model the effects 

That intense winter storms have on the upper ocean thermal structure. 
The largest vertical flux of turbulent kinetic energy (Fig. 3-ID) may be 
expected at this station throughout the season, and makes this an ideal 
location for examining the parameterizations of dissipation postulated 
Beeecis and turner (1967), Kim (1976), and Elsberry, et al. (1976). 
Additionally, it should prove instructive to determine if the major ther- 
mal structure changes may be accounted for with a model that assumes that 
the density structure may be represented by only the temperature struc- 
Ture, and surface buoyancy flux by only the heat flux. 

The study of the processes which control the evolution of the upper 


ocean at OWS V and OWS N offers an opportunity to examine these mechanisms 
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under different conditions than at OWS P. OWS V will be characterized 
by the largest upward surface heat fluxes (Fig. 3-IC) observed at any 
station together with large and variable energy fluxes. Therefore, both 
mechanical mixing and free convection may be expected to play important 
roles in modifying the thermal structure. 

OWS N, on the other hand, is an excellent location for testing the 
model's performance under rather weak and steady forcing. Because of 
the characteristic ocean structure in the subtropics, neglecting salinity 
should not introduce any severe limitations on model performance at 


either OWS V or OWS N (see Dorman, 1974). 
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Reson eeOr Ine UmreRsOOrAN) TOs TRONG 
ATMOSPHERIC FORCING: OBSERVATIONS 
AND SIMULATIONS 
A. |NTRODUCTION 

In this chapter we investigate the oceanic thermal response asso- 
ciated with strong autumn and Satin winter atmospheric forcing events 
at the three ocean weather stations. Denman and Miyake (1973) presented 
data from OWS P (13-24 June 1970) which illustrated the behavior of the 
upper ocean during the passage of several synoptic-scale weather systems. 
Further, they were able to simulate the major features of the mixed layer 
response to these summer storms using a numerical version of the Kraus- 
Turner model (1967). However, the response of the upper ocean to the 
strong fall and winter storm events has not been adequately investigated 
and has never been successfully simulated. 

The atmospheric forcing associated with the strong events that occur 
during the fall and early winter seasons is typically much more energetic 
than the forcing observed during The strong summer events. In addition 
to large downward fluxes of mechanical energy, the fall and winter events 
are frequently characterized by large upward turbulent fluxes of latent 
and sensible heat. Tere esemecheatcal mixing and convection play an 
increasingly important interactive role in the mixed layer evolution 
during the autumn and early winter seasons. 

The principal objectives of this chapter are to examine data sets 
representing air-sea interactions at OWS P, OWS N, and OWS V to: 

|. show that significant upper ocean thermal structure modi fica- 

tions occur during periods of strong fall and early winter forc- 


ing at the three ocean weather stations; 


aa 


2. demonstrate that a modified version of the Kraus-Turner mode! 
(parameterized to account for dissipation enhancement and non- 
penetrative convection) is capable of simulating a large per- 
centage of these changes; 

3. establish the relative importance of the mechanical and convec- 
tive processes that characterize the strong fall! and early 
winter forcing events at the three ocean weather stations. 

In this research the three one-dimensional, bulk meaeis (specifical- 
meeieecr|, and KIM), discussed in Chapter I|, were evaluated at The 
Three ocean weather stations. Successful simulation of the response of 
The upper ocean depends upon adequate parameterization of the principal 
physical processes which govern the response. The relative capabilities 
of the three models in simulating the mixed layer evolutions will! be 
discussed in terms of the importance of properly parameterizing dissipa- 


Tion enhancement and non-penetrative convection. 


Seer ARAMETER! ZATION OF THE MODELS 

A modification of the algorithm presented by Thompson (1976) was de- 
veloped to parameterize the KT, EFT, and KIM models in a consistent 
manner. The scheme was designed to calculate the potential energy changes 
resulting from the vertical fluxes of heat and turbulent Kinetic energy 
That occur in the upper ocean. The sources and sinks of turbulent kine- 
Tic energy, specified in (2-40), are parameterized according to the KT, 
EFT, and KIM models. The modifications of the thermal structure by 
These vertical heat and energy fluxes are accomplished in a manner de- 
scribed in detail in Appendix B. 


oT 


The time rate of change of the mixed layer temperature (<=) and 


depth cn) , specified by (2-21) and (2-33), are specified by the 


3) 





vertical redistribution of heat by the turbulent energy fluxes. In 
Thompson's algorithm the temperature profile is stored in N equally- 
spaced grid intervals [-(n-|) Az, -nAZ], where n=1,2,...,N and AZ = 
2.5m. At each time step (At = | hr) the mixed-layer temperature is 
taken to be the temperature in the first interval. The scheme does not 
eoieitly calculate the mixed-layer depth and it is simply defined Cin 
the mode! and the BT data) as the depth af which the temperature is 
0.2°C less than the mixed layer temperature. 

The basic algorithm was evaluated by Thompson (1976) for the case 
of mechanical mixing with no radiation, and for surface cooling with no 
mechanical forcing. Numerical results from this scheme were equivalent 
to the analytical! calculations for these cases. A desirable property 
of Thompson's algorithm is that it conserves heat and potential energy 
at each time step. This property is important because it allows us To 
isolate and compare the relative importance of mechanically generated 
Turbulent kinetic energy and free convection in the modification of the 


upper ocean thermal structure. 


C. INPUT DATA 

The atmospheric forcing ee Oc, Qe Qn Q,? was estimated from the 
routine meteorological observations using the bulk formulas described in 
Chapter ||| and Appendix A. The turbulent energy fluxes oe, Qo Q,) 
were calculated every three hours and linearly interpolated to hourly 
values corresponding to the integration time-step used in the models. 
The calculated daily value of effective long-wave radiation (Q,) was 
distributed uniformly over the 24 hours. The value of net short wave 
radiation (Q.), however, varied as a function of the.solar altitude. 
To approximate this effect, the estimated daily value was distributed 


Sinusoidally from sunrise to sunset. 
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Since the forcing was calculated from the observed surface and near- 
surface parameters, there was no feedback possible between the ocean and 
the atmosphere. The models were purposely run in this uncoupled state 
to insure that each model received the same forcing. The different model 
responses may, therefore, be compared in terms of the differences in the 
internal physical mechanisms specified by each model. Additionally, since 
any eenors In the forcing are introduced in each model, they should not 
adversely affect the relative comparison of the different model results. 

The models were initialized using mechanical BT data from the NODC 
historical file at the three ocean weather stations. In choosing data 
Sets for testing the performance of the models, the lack of BT data was 
the most limiting factor. During many strong forcing events the BT data 
were completely missing. a eee data sets could only be initialized 
using observations which were available before the events and the results 
of the models validated using data obtained a few days after the events 
had passed. After examining many data sets, it was observed that the 
morning BT (within 2 hrs of 0800) provided a reliable indicator of the 
Trend in the mixed-layer depth and temperature during the strong fall and 
winter events. The individual morning BT observations were examined and 
a number of simple gross error checks were made to insure that the obser- 
vations were realistic. The trend in the sea-surface temperature in the 
BT observations was compared with the trend in the bucket temperature 
reported in the marine deck. If these two trends were not comparable, the 
data set was rejected. Additionally, the trend in the mixed-layer depth 
was examined to insure that it was compatible with the character of the 
computed atmospheric forcing. Finally, to guard against any single BT 
Observation seriously biasing the results (especially the initial obser- 


vation), the profiles used in the model were smooth by hanning. This 
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morning sampling rate prohibits ener of model performance during 
time scales of less than one day. Therefore the analysis will focus on 
the relative capabilities with which the models simulate the general 
trend in the mixed layer changes on time scales of a few days to a few 
weeks. 

Some of the variation in the mixed-layer depth in the individual BT 
observations can be attributed to internal waves with periods ranging 
from the Brunt-Vaisala period (1-10 min in the upper ocean) to the 
local inertial period (between 15.7 hrs at OWS P to 24 hrs at OWS N). 

The atmospheric events investigated in this study were primarily selected 
To include periods with impulsive increases in the wind speed. Inertial 
waves, investigated by Pollard (1970), and Pollard and Millard (1970), 
which are generated by these moving atmospheric disturbances probably 
account for some of the differences in mipecaienee depth between the 

model and data. However, since these internal oscillations have relative 
small effects on the sea-surface temperature, they are easily distinguish- 
able from the changes over several days that are due to the turbulent 
Processes parameterized in the models. 

eemota! Of 49 data sets (20 at OWS P, |I6 at OWS N and [3 at OWS V) 

were finally accepted and modeled using the KT, KIM, and EFT models. 
The scale depth, Z, for the EFT model (see Eq. 2-53) was calibrated to 
six data sets (two from each ocean weather station). A value of Z = 50m 
gave the best fit to the mixed layer depth, in these calibration experi- 
ments, and was used jn all further applications of the EFT model. It 
should be noted, however, that none of the data sets chosen for presenta- 
Tion were among the six calibration sets. 

In this chapter a total of five data sets will be presented (three 


from OWS P, and one each from OWS V and OWS N). They are illustrative of 
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the different atmospheric forcing and oceanic response characteristics 
that were observed at the three ocean weather stations during this study. 
The first four data sets were selected during periods when the mixed 
layer response appeared to be largely one-dimensional. That is to say, 
the local heat budget was maintained and a one-dimensional mixed-layer 
model, adequately parameterized to account for dissipation enhancement 
and partial ly-penetrative convection, was capable of simulating the mixed 
layer response. The relative performance of the three models will be 
examined using these data sets to gain insight into the problem of pro- 
perly parameterizing the important vertical processes. The final data 
Set was selected to illustrate the capability of the one-dimensional 
model during periods when non-local processes an advection) play a 
Significant role in the evolution of the thermal structure. 
De CHARACTERISTICS OF THE ATMOSPHERIC FORCING AND 
BeenNIC RESPONSE AT OWS P 
In this section three data sets from OWS P are examined to illustrate 
the characteristics of the upper ocean thermal response under the fol low- 
ing set of atmospheric and oceanic conditions: 
|. abnormally strong, impulsive mechanical forcing events occurring 
relatively early in the fall when the mixed layer is shallow; 
2. steady atmospheric forcing (equal to the climatological average) 
occurring in the late fall when the mixed layer is deep; 
4. alternate periods of strong and weak mechanical forcing, in the 
early fall, while the ocean is being heated. 
The first example was selected from the period 7-19 October 1954. 
It is representative of the large upper ocean thermal structure modifica- 


Tions that were observed during periods of strong, impulsive mechanical 
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forcing events occurring early in the cooling season (10 of the 20 data 
sets at OWS P had similar forcing). Figures 4=! to 4=2A present the 
relationships between the atmospheric forcing, the observed mixed layer 
response, and the relative performance of the KT, EFT, and KIM models 
in simulating the response. 

In Fig. 4-IA the atmospheric forcing is represented as daily averages 
of the surface stress (TL), the effective insolation (Q,)> the sum of the 
latent and sensible heat (Q7), and the net heat exchange at the sea sur- 
face (Q,)- In this chapter a positive heat flux represents an upward 
flux. To illustrate the relative strength of the forcing during this 
experiment it was determined that 70% of the turbulent kinetic energy and 
81% of the net surface heat fluxes for the month of October 1954 were ex- 
changed during the period 7-19 October (42% of the month). Furthermore, 
comparing the forcing received during October 1954 with the long-term 
mean forcing for October (see Fig. 3=-I) showed that the turbulent kinetic 
energy flux was 120% larger than normal while the net surface heat flux 
was about average during October 1954, 

The surface stress is characterized by three distinct peaks (events) 
centered on II, 14, and |7 October 1954. The net surface cooling (Q-) 
has only one peak on ||! October 1954 which does not occur with the largest 
peak in the stress. The magnitude of Or as specified by (3-2) and (3-3), 
depends to a large degree upon the magnitude of the air-sea temperature 
and vapor pressure differences, in addition to the wind speed. These 
air-sea differences were smaller during the second event than during the 
iest and third events. 

The response of the mixed layer during this period is depicted in 


Fig. 4-1B, along with the response of the KT, EFT, and KIM models. The 
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Figure 4-|, Atmospheric forcing and observed and predicted mixed 
layer response at OWS P during the period 7-19 October 
[954, 
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long-term mean trend (CLIM) is displayed as a basis for comparing the 
mixed layer response. During this experiment the observed mixed-layer 
depth and temperature (30-52 m, 13.9-10.8°C) were abnormally large com- 
pared with CLIM (30-59 m, 13,9-12.6°C). It should also be observed that 
the mixed layer response predicted by the EFT mode! compared favorably 
with the general trend in the data (30-55 m, 13.9-10.9°C). The KIM 
model (30-73 m, 13.9-9.9°C) and KT mode! (30-81 m, 13.9-9.6°C) predicted 
excessive deepening and cooling. 

Figure 4-2 depicts the capability of the EFT model in predicting 
The thermal structure modifications during the period. The agreement 
between the data and model temperature profile is generally quite good 
with the exception of the abnormally large temperature gradient predicted 
below the mixed layer. In an attempt to correct this undesirable 
Property, vertical diffusion was added to the models. The temperature 


change due to diffusion was specified according Ta: 


ot eZ ell 


where A is the vertical diffusion coefficient. A constant value of 
a, amy. an ae as used by Haney and Davies (1976), gave acceptable 
results in the data sets modeled in this study. The diffusion tendency 
was calculated at each time step with the Euler scheme and added to the 
profile after the turbulent mixing processes were calculated. A zero- 
flux condition was assumed at the top and bottom of the profile. The 
resulting accumulation of heat at the bottom was not significant for the 
integrations performed in this study (less than a month in all cases). 


Figure 4-2A shows the results of the EFT model including diffusion. 


it is clear that the mode! with diffusion more realistically represents 
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the structure of the upper thermocline during This experiment. Addi- 
tionally, the diffusion results in a predicted mixed layer which was 
0.1°C cooler and 2 m shallower because of the additional downward heat 
flux. It was concluded that vertical diffusion was an effective means 
for preventing abnormally large temperature gradients, and hereafter 
was included as an integral part of all the models. 

As indicated by (2-21) the temperature of the mixed layer is modi- 
fied by the net heat flux at The oe interface (QO) and the downward 
entrainment heat flux at the base of the deepening mixed layer (w'T'(h) 
= “ho AT). Figure 4-3 depicts the relative magnitude of Q. calculated 
from the marine observations and w'T'(h) calculated from the three models 
(for example BASE/EFT). Also included in this figure is an estimate of 
the entrainment heat flux from the data (BASE/OBS).°  BASE/OBS is the 


change 


before mixing 


~7 =e = eee hoe mixing 


Sea content represented by area A in the schematic, less the contri- 
bution of the surface heat flux (Qo). For instance if Q. > 0 4 the 
contribution during the deepening would be to increase A and therefore 
oh must be subtracted. If 2. <2 Os 00 is added to the heat content in 
area A. If the local heat balance is maintained during this evolution, 
Then this calculation will be representative of the entrainment heat 


flux. Both measurement errors and non-local processes contribute to 


errors in this estimate of entrainment heat flux. 


GZ 





Figure 4-4A indicates that the cumulative entrainment heat flux cal- 
culated by the EFT model (-8350 ly) was in good agreement with BASE/OBS 
(-9200 ly). The KT (-!16,850 ly) and KIM (-15,400 ly) models calculated 
unrealistically large entrainment heat fluxes because of the excessive 
mixed-layer depth predicted by these models. It should be noted, however, 
that the entrainment heat flux, calculated by the EFT model, is much 
larger in magnitude than the surface heat flux (+2400 !y) during this 
period. Therefore, according to (2-21), the principal mechanism account- 
ing for the large sea-surface temperature change observed in the data 
was the downward flux of heat occurring during entrainment. The impor- 
tant implication is that during these large mechanically forced events, 
an accurate prediction of the sea-surface temperature depends, to a 
large extent, on the accurate specification of the entrainment heat flux. 
This iS synonymous, in the models, with an accurate determination of the 
potential energy changes That result from mechanical mixing and convec- 
tion (see Eq. 2-49). 

Figure 4-3B depicts the changes in The potential energy as calculated 
by the EFT model. It should be opserved in Figs. 4=-IB and 4-5A that the 
curves representing S and w'T'(h) are mirror images of the changes 
in potential energy due to mechanical mixing represented in Fig. 4-38. 
weececiantticant to note, however, that in Fig. 4-IB the curve depicting 
e 2 is also similar in shape to a pew n(n wane as (The time rate 
of change of total potential! energy). The similarity in the shape of 
these curves is characteristic of al! data sets in which the vertical 
heat fluxes at the base of the layer dominated the heat budget of the 
mixed layer. 


The differences in the performance of the EFT, KT and KIM models may 


be examined in terms of the potential energy budgets calculated by each 
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Figure 4-3. (A) Relative importance of the vertical heat fluxes and 
(B) Potential energy modifications calculated by the 
EFT mode! during the period 7-19 October 1954, 
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model. Table 4-1 presents a comparison between the total surface heat 

fl Ux (20) and the total entrainment heat flux (£ w'T'(h)) estimated 

from the data and from the model results. The final three columns repre- 
sent the total potential energy change (ZAPE), the total potential energy 
increase due to mechanical mixing (LAPE ), and the total reduction in 
potential energy due to convection (LAPE_). These three quantities are 


related as follows: 


TAPE = Se eee ete) ee a aes 


where r is the fraction of convectively-generated turbulent kinetic 
energy available for entrainment and is taken as 0.15. 

Comparing ZAPE for the three models it may be observed that in 
The KT and KIM models, considerably more turbulent kinetic energy is 
available for entrainment than in the EFT model. I+ should also be noted, 
comparing ZAPE with LAPE » that convection was not as important as 
mechanical mixing during this period. An examination of LAPE Shows 
That the largest reduction in potential energy is calculated by the KT 
model because the mixed layer is deepest and therefore convection occurs 
over a deeper depth (see Eq. 2-51). The net result is that ZAPE is 
largest in the KT model and smallest in the EFT model. 

It may easily be Romer erated that the abnormally large mixed-layer 
depth and temperature changes predicted by the KT and KIM models are due 
To the large amount of turbulent kinetic energy available for entrainment. 


One may express (2-33) as 
ae) (4-3) 


where Ab = OP gaAT : 
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Thus the response to mechanical mixing, (ae) » Is the ratio of The tTur- 
bulent kinetic energy available for entrainment (G,-D,-S, = ZAPE) TO 

the amount of work necessary to entrain the denser fluid from below The 
mixed layer (hAb). An examination of the models indicated that hAb in- 
creased throughout the period |2-I9 October, and that if was largest in 

the KT model and smallest in the EFT model. Therefore, the excessive 
mixed layer depth response in the KT and KIM models relative to The EFT 
Ree and the data, was primarily ite to the numerator of (4-3). The 

lack of dissipation enhancement in the KT model, and The linear parameteri- 
Zation in the KIM model, are inadequate during this experiment. However, 


oT 
The data (a0 eric in Fig. 4-18) suggest that the exponential para- 


oT 
meterization of the EFT model was reasonable during these events. Compar- 
ing ZAPE between the KT and EFT models shows that the EFT model calcu- 
lates a dissipation rate equal to 60% of the surface production rate. In 
the KIM model, however, dissipation and storage were only 35% of the 
Surface production. 

The mixed-layer depths and temperatures, predicted by The EFT model, 
were quite good, and & w'T'(h) estimated from the data was in close 
agreement with £ w'T'(h) calculated in the EFT model. This suggests - 
That The potential energy changes calculated by the EFT mode! should be 
reasonable compared with the potential changes occurring in the upper 
ocean (which are not easily calculated over short time periods). 

If is important to recognize the important contribution of the en- 
Trainment heat flux during these large mechanically forced events. The 
sea-surface temperature is reduced during the fall and winter not only 
by surface cooling but also by the vertical re-distribution of the heat 


by mechanical mixing. The implication to modeling the upper ocean response 
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is that an accurate prediction of the changes in the mixed-layer depth 
is essential for an accurate prediction of the sea-surface temperature 
response. 

The second case to be examined was chosen during the period 22 
November-!2 December 1957. It is representative of the upper ocean 
thermal response that occurs later in the season when, relative to the 
Previous example, the mechanical forcing is weaker and less impulsive. 
In contrast to the previous case, however, the net surface cooling that 
occurred during this period was larger and, as will be feoneanonee 
plays a more important role in the evolution of the mixed layer. The 
Surface stress, depicted in Fig. 4-4A, is rather steady during this 
Period, with the exception of two weak events centered on 26 November 
and 6 December 1957. The net surface cooling is also steady with no 
large peaks. Comparisons between ae and Q- with the climatological 
trend (Fig. 3-1) showed that the transfer of turbulent kinetic energy 
during this period was about average (104% of the normal value). How- 
ever, more surface cooling occurred during this period than the average 
(120 ly/day more than climatology). The mechanical forcing in this 
data set is typical of 9 of -the 20 data sets investigated at OWS P. 

If may be observed from DATA in Fig. 4-48 that the response of the 
upper ocean thermal! structure (58-78 m, 10.1-8.7°C) was very close to 
the long-term mean trend (58-77 m, 10.1-9.0°C). It should also be 
observed that this steady response was closely simulated by the EFT 
model) (56-80 m, 10.1-8.6°C) while the KT (58-118 m, 10.1-8.0°C) and KIM 
(58-104 m, 10.1-8.0°C) models again predicted excessive deepening and 
cooling rates. The structure modifications predicted by the EFT model, 


shown in Fig. 4-5, also agree very closely with the data. The model not 
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only predicted approximately the correct mixed-layer temperature and 
depth, but accurately simulates the reduction in the temperature gra- 
dient (AT) Immediately below the mixed layer. As shown by Gill and 
Turner (1976) the reduction of AT during these deepening events is pri- 
marily caused by non=-penetrative convection. Equation (4-3) may be 


examined to verify that a reduction in AT results in the reduction of 


dh) 


37m ° The reduedi onalnmecdaee Ss 


hAb and, therefore, an increase in ( 
the principal mechanism by which convection interacts with mechanical 
mixing to permit further mechanical epenthe of the mixed layer. 

In Fig. 4-6A and Table 4-2 it should be observed that during this 
period the total net surface heat flux (+5497 ly) was considerably 
larger than the entrainment heat flux estimated from the data (-2972 ly), 
or calculated by the EFT model (-348! ly). Note also that the entrain- 
ment heat flux calculated by the EFT model is again in good agreement 
With The estimate from the data. The KT (-10240 ly) and KIM (-8132) 
models again calculated abnormally l!arge entrainment fluxes (because of 
inadequate dissipation) eee a consequence overestimated a 

The potential energy change observed in Fig. 4-68 and Table 4-2, 
calculated by the EFT model, is very different than in the previous data 
set (see Fig. 4-2B and Table 4-1). With the exception of the initial 
increase during the first few days, the general trend is for the poten- 
tial energy to decrease during the period. As a result the characteris- 
Tic shapes of the curves representing ae a and 2 are not 
Similar as in the previous example. This property is characteristic of 
The EFT model when convection plays a dominant role in the evolution of 


the upper ocean. The potential energy changes calculated by the KT and 


KIM models may be examined in Table 4-2 to see that a total increase in 
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potential energy resulted in both models. However 4 comparison of 
oT 


gh = , and £ w'T'(-h) between the models and the data suggests 
nat oo calculated by the KT and KIM model is unreasonable. These 
comparisons further indicate that ss calculated by the EFT model 
should be representative of as occurring in the upper ocean during 
this period. 


The excessive deepening and cooling ratio predicted by the KIM and 
KT models may again be related to an inadequate parameterization of 
dissipation enhancement. In the EFT model the dissipation rate was 75% 
of the surface production rate while in the KIM model the dissipation 
and storage rate was only 55%. This is to be compared with the 60%/ 
35% calculated in the EFT/KIM model in the previous example. The in- 
creased dissipation during this second case is the result of a deeper 
mixed~layer in this example (58-78 m) relative to the previous case 
(30-52 m), and the dissipation parameterizations defined by (2-53) and 
(2-53A). As with the previous example the large LAPE calculated by 


The KT and KIM models result in a serious overestimate of au and 
aE 
ot ° 

I+ should be observed in Tables 4-| and 4-2 that for the EFT model 


LAPE was smaller in the second example than in the first. However, 


it may be observed in Figs. 4-IB and 4-4B that a is comparable in 
aT 
both cases but —— was much smaller in the second case. One can see 


from Eq. (4-3) that since G,-D, is smaller in the second case, hAb must 
also be smaller for o to be comparable. Since h was also larger 


in the second case, the gradient at the base of the mixed layer must be 
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smaller, which may be confirmed in Figs. 4-3 and 4-5. If we simply ex- 


press (2-20) as 


oT LO + Dae h } 
ae id 


ae h (4-4) 


then the smaller decrease in the sea-surface temperature during the 
Second case, relative to The first, is understandable. Comparing BI 
and LZ w'T'(-h) in Tables 4-1 and 4-2 shows that the sum of the net 
surface and entrainment heat fluxes were smaller in the second case and, 


as previously mentioned, the mixed-layer depth was larger. Equation 


oT 
(4-4) would therefore explain the smaller == in the second case. 
Equations (4-3) and (4-4) indicate that, for comparable forcing, a 
oT 
larger mixed layer response (ae and =) will take place when the 


layer is shallow than when it is deep. The important implication is 
That the large storms which take place early in the cooling season wil! 
have a much larger effect on the mixed-layer depth and temperature than 
Those occurring relatively later. 

The final data set to be presented in this section was selected for 
The period I-15 September 1966. I+ is illustrative of conditions, in 
early September, when the mechanical forcing is characterized by alter- 
nate periods of relatively strong and weak forcing, and when net surface 
heating is occurring. Under these conditions the mixed layer frequently 
deepens during the strong forcing and retreats during the weak forcing 
periods (see Chapter I1!, Eq. (2-43) for a discussion of this process). 
It will be demonstrated that a proper parameterization of dissipation 


is essential to simulate this type of ocean response. 


1) 





An examination of Fig. 4-/A shows that the mechanical forcing is 
characterized by three events, centered on September 4, 8, and 12, and 
the beginning of a fourth event on September 14. Further it should be 
observed from Fig. 4-7A and 4-9A that net surface heating takes place 
Throughout most of this period, while the mixed-layer temperature is 
decreasing (Fig. 4-7B). Because net surface heating is occurring while 


dT 
——< 0, (4-4) shows that the decrease in the sea-surface temperature 


dt 

is entirely due to the entrainment heat flux. Therefore during periods 
characterized by net surface heating, the sea-surface temperature change 
is predictable only if the entrainment process is properly parameterized. 

In response to this alternately weak and strong mechanical forcing, 
the predicted mixed layer depth alternately retreats and deepens (Fig. 
Seoveeerigures 4-/B and 4-8 illustrate that the EFT model simulates the 
evolution of the mixed layer quite well during this period. As expected, 
the KT model consistently predicts a deeper, cooler mixed layer than the 
KIM and EFT models. However, during this period, contrary to the previous 
examples, the KIM model predicts a warmer and shallower mixed layer than 
the EFT model. The exception to this is during the periods following 
the first and second events. This Suggests that less turbulent kinetic 
energy is available for entrainment in the KIM model than in the EFT model, 
except during the peaks in the forcing. 

Table 4-3 may be examined to show that the three different predicted 
mixed layer responses are due to the different assumptions regarding 
dissipation enhancement. Comparing LAPE and LAPE it may be observed 
That convection played a relatively minor role during this period. The 
interesting observation is that LAPE is ema in the EFT model than 


in the KIM model (as the observations from Fig. 4-7B suggest). This is 
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a departure from the previous data sets and may be understood by an 
examination of the dissipation enhancement parameterization in each ° 
model. 


We may express (2-53) and (2-53a) in the following form 


Z = 0 for Ki 


= 50m for EFT (4-5) 
Dox 
rea Doh 
= “|e acs (4-6) 
Ae: Wy 


Equations (4-5) and (4-6) represent the ratio of turbulent kinetic 
energy (G.-D,) available for mixing to The surface production accord- 
ing to the KT model (D.Wy) The ratios in (4-5) and (4-6) are 
plotted in Fig. 4-I0 as a function of various wind speeds and layer 
depths. They may be interpreted as the percentage of mechanical pro- 
duction of turbulent energy (relative to surface production of the KT 
model) available for entrainment in the various models. In the KT 
model dissipation enhancement is not parameterized and G,-D, is always 
equal to OW Seen ine eh) model) Gl=by wais@oeual io OW only 
for the trivial case of h=0. In general G,-D, < OW according 
To (4-5) and therefore the KT model always predicts a deeper, cooler 
mixed layer than the EFT model. However, because the surface produc- 
tion in the KIM model was parameterized according to the Kato and 
Phillips (1969) experimental results, it is 25% larger than in the KT 
and EFT models. Dissipation enhancement is parameterized in the KIM 


model as a linear function of depth and a constant background 
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Figure 4-10. Vatues G,-D, normalized by op ‘ee dS a Tuner ionncek 
depth and wind speed for the EFT, KT, and KIM models. 
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dissipation (see Eq. 2-53a). In Fig. 4-10 it should be observed that 
there are combinations of wind speeds and mixed layer depth for which 
the amount of kinetic energy available for entrainment predicted by the 
KIM model may be greater or smaller than either the KT or EFT models. 
During the events characterized by strong mechanical forcing, the 
response of the KIM model is very similar to the KT model and excessive 
mixed-layer depths are predicted. However, during periods of weak 
mechanical forcing, the mixed-layer depths predicted by the KIM model 
are too shallow and consequently the mixed-layer temperatures are too 
warm. These model results further suggest that the exponential para- 
meterization of G,-D, , in the EFT model, may be calibrated for a 
wider range of wind speeds and layer depths than the linear parameteri- 


Zation in the KIM model. 


E. IMPORTANCE OF NON-PENETRATIVE CONVECTION 

In the previous three examples the percentage of convectively- 
generated turbulent kinetic energy available for entrainment was set 
equal to 15% (r = 0.15 in Eq. 2-50) following Gill and Turner (1976). 
In This section it will be demonstrated that non-penetrative convection 
(cooling of the mixed layer without deepening) should be an integral 
feeneon any turbulent bulk model. Furthermore, this property is 
absolutely essential during the fall and winter cooling seasons to pre- 
vent the models from predicting excessive deepening and cooling rates. 
Gill and Turner reported that the potential energy of the upper 250 m 
at nine ocean weather stations in the Atlantic reaches a maximum by 
early October and decreases during the October to March period. They 


were successful in simulating these potential energy changes with a 


So 





* 





modified version of the KT model with r = 0.15. The experiments of 
Deardorff, et al. (1969) and Farmer (1975) suggest that r is of the 
order 0.0I1-0.15 but there were uncertainties in their estimates. How- 
ever, the essential feature that must be simulated by the models is a 
decrease in the total potential energy during periods when the upward 
surface heat flux is larger than the entrainment heat flux. 

The data set chosen to illustrate this point is from OWS V during 
the period 7-19 October 1963. It was selected because it is typical of 
the response of the upper ocean at OWS V during periods when the net | 
surface heat flux is positive and larger than the entrainment heat flux 
(10 of the 13 data sets examined at OWS V are in this category). In 
This experiment, only the EFT model was used and the fraction of convec- 
Tively-generated turbulent kinetic energy available for entrainment was 
varied between 15% and 1|00%. 

Figures 4-|| to 4-13 depict the calculated forcing and the relative 
capability of the EFT model with varying amounts 5 non-penetrative 
convection. The atmospheric forcing (Fig. 4-I!A) is characterized by a 
peak in the stress and surface cooling (Q7) centered on 14 October. Com- 
peared with the forcing in the first example presented (Fig. 4-IA), the 
Turbulent kinetic energy flux was nearly five times smaller, but the 
Meteourface heat flux was larger by 550 ly. In the next chapter it wil! 
be demonstrated that large upward heat fluxes are common at OWS V because 
of the extremely large air-sea temperature and vapor pressure differences 
during the fall and early winter seasons. 

In Fig. 4-11B it should be observed that a significant mixed layer 
response occurred (42-60 m, 24.3-23.5°C). It should also be noted that 


The EFT model (r = 0.15) simulates this response with excellent agreement 
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Figure 4-13, Same as Figure 4-3 except for OWS V during the 
period 7-19 October 1963. 
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with the data (42-61 m, 24,.3-23.4°C). Note also in Fig. 4-12 the good 
agreement between the EFT model (r = 0.15) profiles and the data pro- 
files. However, the EFT model with r = 0.5 (42-68 m, 24.3-23.2°C) and 
r= 1.0 (42-74 m, 24.3-23.0°C) predicted deepening and cooling rates 
that were too large. 

In Fig. 4-I3A and Table 4-4 it should be observed that the net surface 
heat flux (+3000) was larger than the entrainment heat flux estimated 
from the data (-1870). As previously mentioned, this should result in 
a total decrease in the potential energy and Fig. 4-138 and Table 4-4 
show that the EFT model with r = 0.15 predicts this response. An examina- 
Tion of the potential energy modifications by the three variations of the 
EFT model (Table 4-4) illustrates the importance of properly parameteriz-~ 
ing non~penetrative convection. It is important to observe that for 
r= 0.15 (15% of the convectively-generated turbulent kinetic energy 


available for entrainment) and r = 0.5 the total potential energy (2ZAPE) 


is decreased in the model, while for r = 1.0 it is increased. In fact, 
as discussed in Chapter || (see 2-51), a model with full penetrative 
convection (r = 1.0) must increase the potential energy. In accordance 


with the findings of Gill and Turner (1976) a model with r = 1.0 will 
not be useful for seasonal integrations. The present results from the 
EFT model with r = 1.0 indicates that, in certain circumstances, it wil 
also be wrong for very short integrations during the cooling season. 
The results from the EFT model with r = 0.15 suggest this model (with 
an exponential dissipation parameterization) reasonably simulates this 
Process. 

An additional reason for presenting this data set is because its 


time frame (7-19 October) is the same as the first example (Fig. 4-1) 
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from OWS P. Since these two examples are typical of the atmospheric 
forcing events and associated oceanic response at these two stations, 
they may be compared to illustrate the essential differences in the 
physical mechanisms that modify the upper ocean at these two stations. 
Comparing LAPE (for EFT, r = 0.15) in Tables 4-| and 4-4 it is ob- 
served that nearly four Times as much turbulent kinetic energy was 
available for entrainment at OWS P as af OWS V. It should be further 


noted that the deepening rates coi were comparable but resulted in 


ot 
a significantly larger decrease in the mixed-layer temperature at OWS P 
oT 
(compare 37 and Sa in Figs. 4-IB and 4-11B). Considering The 


ratio defined by (4-1) hAb must be significantly larger at OWS P. 
Since the mixed layer is deeper at OWS V during these data sets the 
below layer gradient (represented in the mode! as AT) at OWS P must 
account for hAb being larger at OWS P. It may be observed in Figs. 
4-2A and 4-|2 that the below layer gradient at OWS P is typically 

much larger than at OWS V. Therefore Significantly more work is re- 
quired to deepen the layer at OWS P. However, for comparable deepening 
rates as in these two cases, the resulting entrainment heat flux wil 
be much larger at OWS P than at OWS V (compare 2 w'T'(-h) in Tables 
4—-| and 4-4 during these two cases). The larger entrainment flux at 
OWS P is the principal reason for the larger sea-surtace temperature 
change, since the net surface heat fluxes (2Q_) are very similar during 
These two cases. Since these two examples are typical of the response 
of the mixed layer at these two stations during the fall and winter 
storm events it is suggested that the entrainment heat flux plays the 
most dominate role at OWS P. While the entrainment heat flux at OWS V 
is less important, it is still an essential part of the local heat bud- 


get of the mixed layer, and can not be ignored. 
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In this and the previous section the data sets examined illustrated 
the different upper ocean responses that occur under various atmospheric 
forcing and initial oceanic conditions. It should be noted that large 
upper ocean thermal responses may take place when either the mechanical 
forcing or upward turbulent heat fluxes are large. It is also important 
to note once again that both the surface heat flux and the entrainment 
heat flux must be specified in order to predict the sea-surface tempera- 
ture. It should also be noted that during these large storm events, the 
response BA the ocean is primarily one-dimensional because the vertical 
heat fluxes are extremely large during relative short periods. To ade- 
quately simulate these large One-dimensional responses it has been 
demonstrated that an adequate parameterization of dissipation enhance- 
ment and non-penetrative convection is absolutely essential. A combina- 
Tion of r = 0.15 (after Gill and Turner, 1976) and the EFT model para- 
meterization of dissipation enhancement consistently were in better 


agreement with the data than the KT or KIM models. 


Seer ORIANCE OF NON LOCAL EFFECTS 

In this section the performance of the three models will be examined 
during a period when the local heat balance was not maintained. The 
Purpose of this presentation is to demonstrate that the one-dimensional 
model is capable of providing useful temperature structure information 
during periods characterized by large horizontal heat fluxes. The data 
Set chosen to illustrate these points is from OWS N during the period 
9-24 November 1965. Figures 4-14 to 4-16 again depict the nature of 
The forcing, the observed mixed layer response, and the relative perform- 
gace of the modeis. The data set is characterized by two distinct events 


centered on November 15 and 22 (see Fig. 4-14A) and the turbulent kinetic 
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energy and net surface flux exchanged during These events were very 
strong compared with climatology (Total ie for November 1965 was 205% 
of the climatological mean for November and total Q_ was 160% of the 
mean). During the period 9-29 November 1965, 95% of the total monthly 
turbulent kinetic energy and 80% of the net surface heat flux were 
exchanged. 

It-should be observed in Fig. 4-14B and 4-15 that the EFT model 
accurately simulates the evolution of the mixed layer through November 
19 and suggests that the local heat balance was maintained during This 
period. The estimate of cumulative entrainment heat fiux from The data 
(Fig. 4-I15A) on November I9 also supports this suggestion. However after 
November 19-20 the observed mixed-layer temperature decreased signifi- 
cantly, relative to all the models, while the mixed-layer depth was 
Simulated quite well by the EFT modei. Furthermore, in Fig. 4-16A if 
may be observed that the entrainment heat flux, estimated from the data 
from November 19-29, appears to be too large. The layer depth, accord- 
ing to the data, only increased 8 m during the period November 19-29. 
This smal | a coupled with the weak temperature gradient immediately 
below the mixed layer (Fig. 4-15), would not support BASE/OBS (Fig. 4- 
ISA) during this period. 

The reason for the difference between the observed and model mixed 
layer temperature during November 20-29 is presumed to be due to a hori- 
zontal intrusion of a cold water mass into the region. This presumption 
iS based upon the observations from Fig. 4-15 that the entire data pro- 
file (even below the mixed layer) is colder relative to the EFT model. 
During this period the ocean weather ship reported its position to within 
1O NM of 30N - I40W and therefore ship drift was not responsible for the 


vemperature changes. Since the temperature changes observed in Fig. 4=\|\9 
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can not be due to vertical mixing, it is suggestive of a horizontal 
heat flux. 

During this period the EFT model was able to simulate the mixed- 
layer depth quite well and the characteristic slope of the temperature 
profile was maintained. If the principal motivation for this forecast 
was to predict the thermal structure modifications, say for predicting 
the resulting changes in the acoustic properties, then the EFI model 
results would be very useful. These model results further suggest that 
if the horizontal heat fluxes were specified, say from an ocean general 
circulation model, then the EFT model would be capable of predicting 
the total response of the upper ocean with adequate accuracy. 
eee LAl IVE PERFORMANCE OF THE EFT, KT, AND KIM 

MODELS AT OWS P, OWS N, AND OWS V 

In this section the relative performance of the EFT, KT, and KIM 
models will be evaluated at the three ocean weather stations. Figures 
4-|7 to 4-19 depict the daily observed mixed-layer depths versus the 
predicted depths for the three models at the three ocean weather sta- 
tions. The diagonal lines are drawn to aid in relating model results 
To the locus of perfect prediction. In Table 4-5 the RMS and mean 
errors for the mixed-layer temperature (MLT) and mixed-layer depth (MLD) 
are presented for the three models at the three stations. A negative 
mean error in MLD indicates that the model, on The average, predicted a 
deeper MLD. A negative MLT indicates that, on the average, the model! 
Predicted MLT was warmer than the observations. The RMS and mean errors 
are based on the daily (morning) observed and predicted MLT and MLD and 
are Dased on 313 values at OWS P, 238 at OWS N, and 208 at OWS V (this 


Includes all 49 cases studied). 
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Figure 4-17, Predicted vs. observed mixed layer depths (MLD) 
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Figure 4-18. Same as 4-!7 except for OWS 
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TABLE 4—5. 


Model 


Model 


Model 


OWS P 
MLT 
MLT Mean 
RMS (°C) Error (°C) 
0.36 +0.16 
0.49 +0. 32 
OL7Z +0.54 
OWS N 
MLT 
MLT Mean 
RMS(2@) Errern(cc) 
0.34 —~0.01 
0235 —~0.02 
0.37 +0.04 
OwS V 
MLT 
MULT Mean 
RMS ( °C) (Errore) 
0265 +0.02 
0.66 +0.05 
0275 +0.18 


MLD 
RMS (m) 


ows 
13.8 
Uoe2 


RMS (m) 


6.7 
Beg 
10.0 


RMS (m) 


ae 
Ee 
14.1 


Comparison of mean and RMS errors of the models 
at the three ocean weather stations. 


Mean 
Error (m) 


oe 
= 136 
=—6 


Exror (m) 


+1.6 
+2.2 
—2.8 


Error (m) 


Ss 
—4.2 
Ie 


* 
RMS and MEAN ERROR based on comparison of daily observed and predicted 


mixed-layer depth (MLD) and mixed-layer temperature (MLT). 


Statistics 


are computed from 313 comparisons at OWS P, 238 at OWS N and 208 at 


OWS V. 
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1+ should be observed in Table 4-5 that the RMS and mean errors for 
MLT and MLD are consistently smallest for the EFT model and largest for 
the KT model. At OWS P, for example, the scatter observed in Fig. 4-1/7 
and the mean MLT and MLD errors in Table 4-5 show that the models are 
all biased toward predicting a mixed layer that is too cool and too 
deep. However the RMS and mean errors, and the magnitude of the scatter 
in Fig. 4-17, demonstrate clearly that the EFT model simulates MLT and 
MLD more accurately at OWS P. During the fall and early winter, OWS P 
is characterized by the largest soeraniiGall forcing observed at the 
Three stations (see Fig. 3-IC). The performances of the EFT model rela- 
Tive to the KT and KIM models suggest that the exponential parameteriza- 
Tion of dissipation enhancement (EFT) was most effective in preventing 
excessive deepening. 

At OWS N the mean errors in MLT and MLD, observed in Table 4-5, and 
The scatter in Fig. 4-18 indicate that the EFT and KT predict MLT and 
MLD, on the average, too warm and too shallow. The KT model is again 
biased toward predicting MLT and MLD too deep and too cool. At OWS N 
considerably less mechanical energy is transferred to the ocean Than at 
OWS P (again see Fig. 3-IC). Therefore, at OWS N, one would expect con- 
vection to play a more dominate role in the mixed-layer evolution and 
The performance of the three models to be more similar. The similarity 
in the scatter in Fig. 4-18 and the RMS and mean errors in Table 4-5 
should be observed to beat Similar at OWS N. An interesting observa- 
TiOn may be made by comparing the relative magnitude of the RMS and mean 
errors at OWS P and OWS N. At OWS P the relationships between the mean 
and RMS MLT and MLD errors suggest that the increasing RMS errors are 


largely due to the increasing bias of the EFT, KT, and KIM models. At 


On 





OWS N, however, the large RMS errors, relative to the mean errors (most 
obvious in MLT), would suggest that the errors are more random at OWS N 
than at OWS P. 

At OWS V the scatter observed in Fig. 4-19 and the RMS and mean 
errors in Table 4-5 indicate That the models are all biased toward pre- 
dicting MLD and MLT too deep and too cool. Again if should be observed 
that the scatter in Fig. 4-19 and the RMS and mean error for MLD are 
smallest for the EFT model and largest for the KT model. As indicated 
previously, mechanical energy and convection are strong at OWS V. The 
increasing RMS and mean MLD errors for the EFT, KT, and KIM may once 
again be related to the parameterization of dissipation enhancement in 
each model. Again it should be observed that the relationship between 
The RMS and mean errors in MLT would suggest that a considerable amount 
of the RMS MLT error is random. 

The data presented in Figs. 4-I/7 to 4-19 and Table 4-5 clearly indi- 
cate that the EFT model is capable of predicting a larger percentage of 
the observed MLT and MLD changes observed in the data sets modeled in 
This study. The exponential parameterization of dissipation enhancement 
coupled with an adequate estimate of non-penetrative convection appears 
To be most reasonable for predicting the changes associated with large 
atmospheric forcing events. 

Table 4-6 presents a comparison of the accuracy with which the EFT 
model predicted the total observed MLT change during the 49 data sets 
examined in this study. The values in Table 4-6 represent the number 
of data sets during which the EFT model predicted the various percentages 
of the total observed MLT change. For example, at OWS P, in 7 out of 
The 20 data sets, greater than 90% of the observed MLT change was pre- 


dicted by the EFT model. I+ may be observed in Table 4-6 that the EFT 


eZ, 





TABLE 4-6. Camparison of MLT predictions by the EFT model at the 
ocean weather stations. 


Percentage of MLT Change Predicted 


Number of 

OWS >90% >75% >50% >25% Data Sets 
P ni AZ 18 18 20 
N 2 i 12 14 16 
V 6 7 10 a 13 


*Values in the table are the number of data sets which predict the 
prescribed percentage of MLT. 
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model showed reasonably good capability in predicting the MLT. For ex- 
ample, in nearly 50% of the data sets modeled, the EFT mode! was capable 
of predicting greater than 75% of the observed MLT changes. Additionally 
Soro Out of 20 cases at OWS P, I2 out of 16 at OWS N, and [0 out of 13 
at OWS V the one-dimensional mode! accounted for greater than 50% of the 
observed MLT change. It should be noted that in this study only a mini- 
mum of calibration was performed with the EFT model. The performance of 
_ the model was still quite good and suggests That greater reliability 
could be obtained with more extensive calibration. 

In addition to providing reasonable estimates of MLT and MLD changes, 
the one-dimensional model provides a means of separating the contribu- 
tions of the vertical heat fluxes to the total heat budget of the mixed 
layer. As may be seen in Fig. 4-I6A (BASE/OBS), The estimation of the 
entrainment heat flux from the data is very unreliable during periods 
when non-local processes are important. In the next section we will ex- 
amine IQ and Zw't'(-h) from the EFT mode! results in an attempt 
To determine the relative importance of these heat fluxes at the three 
weather stations. 

SeeeeeeeAl IVE IMPORTANCE OF THE SURFACE 

AND ENTRAINMENT HEAT FLUXES 

In this section an attempt wil! be made to determine the relative 
importance of the surface and entrainment heat fluxes during the periods 
modeled in this study. In the first four examples presented in this 
chapter, it was demonstrated that during periods when the local heat 
budget was maintained, the EFT model predicted the mixed layer depth and 
Temperature quite accurately. Furthermore estimates of the entrainment 


heat flux from the data were in good agreement with mode! calculations. 
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Since the thermal structure modifications calculated by the model agreed 
closely with the data, it Is assumed the potential energy changes calcu- 
lated in the model are representative of the potential energy changes 
occurring in the ocean. Additionally it will be assumed that during 
meleds wnen non-local processes are important the vertical turbulent 
heat fluxes are still represented accurately by the EFT model. 

Figure 4=20A depicts the relative magnitudes of the total cumulative 
net surface and entrainment heat fluxes calculated by the EFT model. 

In Fig. 4-20B the total potential energy change (ZAPE) is plotted against 
the total upper ocean heat content change (20) calculated by the EFT 
model. Since the model is one-dimensional AH = 0, - the data in 

Figs. 4-20A,B is representative of the total changes calculated by the 
Sedemsin aii but three data sets at the ocean weather stations. One data 
set at OWS P, presented earlier in this chapter, and two at OWS V were 
examples of weak surface heating periods and the entrainment heat flux 
totally dominated the heat budget of the mixed layer. 

These two figures may be Interpreted in light of the four one- 
dimensional data sets previously presented. If ZAPE > O (Fig. 4-20B) 
during these periods, then the mechanical mixing by the wind was the 
dominant vertical turbulent process, and normally this will be indicated 
in Fig. 4-20A by |Z w'TTCh) > [ZQ |. If, however, ZAPE < 0 then 
convection was the dominant vertical turbulent process and |Z w'TT(-h)| 


< [2Q, 





- The relationships depicted in Figs. 4-20 indicate that the 
atmospheric forcing events examined at OWS P were mainly dominated by 
mechanical mixing. In [5 of the 19 cases examined, |Z w'T'Ch)| > }ZQ_| 
and Bae > O . At OWS N and OWS V, however, the large majority of the 


atmospheric forcing events were dominated by convection. In 12 of the 
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Figure 4-20. 


Relative magnitude of the cumulative heat fluxes and 

total potential energy change at the three ocean stations. 
(A) Cumulative net surface heat flux (2Q,) vs. cumulative 
entrainment heat flux (w'T'(-h))(B) Total potential 
energy change (ZAPE) vs. ocean heat content change (AH). 


All points calculated by EFT model for total duration of 
data set. 
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|6 data sets at OWS N [Z w'T"(-h)| < [ZQ_| and ZAPE < 0 , and at OWS V 
10 of Il cases indicated |Z w'T'(-h)| < }zQ_ | and ‘ZAPE < Oi 

It has been demonstrated that the large forcing events, whether domi- 
nated by mechanical forcing or surface cooling, are capable of causing 
extremely large changes in the mixed layer in relatively short periods. 
lt has further been illustrated that the response during these large 
events is largely one-dimensional. Additionally, it has been shown that 
a one-dimensional model, properly parameterized to include dissipation 
enhancement and non-penetrative convection, is capable of predicting a 
large percentage of the observed changes. In the next chapter we wil! 
attempt to determine how significant is the role of the large forcing 


events in the total seasonal evolution of the upper ocean. 
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V. ROLE OF STRONG ATMOSPHERIC FORCING 
EVENTS IN THE SEASONAL EVOLUTION OF 
THE SUPEER OCEAN 
A. INTRODUCTION 
ln the previous chapter it was demonstrated that significant changes 
may take place in the upper ocean thermal structure in response to strong 
atmospheric forcing events. Further it was shown that the one-dimensional! 
processes, when modeled properly, are capable of predicting a large per- 
centage of these observed changes in the three locations studied. In 
This chapter the historical series of surface and near-surface marine ob- 
servations will be examined in a new and rather unique way, and the 
Principal objectives will be to: 

|. determine the significant characteristics of the atmospheric forc- 
ing during the fall and early winter cooling season at the three 
ocean weather stations; 

2. demonstrate that the major features of the upper ocean therma| 
response during the cooling season are explainable in terms of 
One-dimensional processes; 

5. quantify the relative importance of strong atmospheric forcing 


events to the total evolution of the upper ocean at OWS P, 
OWS N, and OWS V. 


Eee HARACTERISTICS OF THE MARINE ATMOSPHERE 

AT THE OCEAN WEATHER STATIONS 

To determine the distribution and variability of the atmospheric force- 
ing at the three ocean weather stations, values representative of wind 
speed (u,), turbulent kinetic energy flux ey and upward turbulent 
heat flux (Q.) were computed from every available three-hourly record. 
These values were grouped into equal (32) class intervals (ranked in order 


Of increasing values) and the resulting frequency distributions are 
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Semered in Figs. 5-| and 5-2, Additionally Table 5-!| lists the signi- 
ficant statistical quantities which characterize these distributions. 
In all future discussions a negative heat flux will represent a heat 
loss by the ocean. 

The non-gaussian nature of these distributions is evidenced by their 
Characteristic shape and the values of skewness and kurtosis presented 
in Table 5-I. It will be demonstrated that the large, relatively rare 
values (reflected by the long tails in the distributions) account for a 
considerable amount of the total energy exchange. The means of the uy, 
and ee distributions show that the strongest winds and largest turbulent 
Kinetic energy exchanges occur at OWS P while the weakest mechanical 
interactions are observed at OWS N. From the distribution of Q. it is 
Btesrved that the largest turbulent heat fluxes Taeonoieee at OWS V and 
the smallest at OWS P. Additionally, Table 5-| shows that the mean of 
Siemeeeaistributton, at all stations, is approximately ABE lie its stan- 
dard deviation. This relationship is also vances Q. at OWS N, while 
eemewo F and OWS V the distribution of 0. has more relative variance 
Than ux. This implies that the variability of the turbulent heat fluxes 
is more closely coupled to the variability of the wind at OWS N than at 
the other stations. It will be demonstrated that the increased variance 
in 0. eee OWS P and OWS V is the result of a farger variability itn the 
air-sea temperature and vapor pressure differences at these stations. 
These distributions indicate that the characteristics of the marine atmos- 
Pheric forcing are similar at OWS P and OWS V but quite different at 
OWS N. The similarities and differences in these characteristics may be 
explained by considering the principal air-sea parameters in terms of 


the geographical locations of the three ocean weather stations. 
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Figure 5-!. Histograms of wind speed (uy) and turbulent kinetic 
energy flux (uy?) at OWS P, OWS N, and OWS V. 
Vertical dashed lines represent the mean of The 
distribution. 
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Figure 5-2. Histograms of upward turbulent heat flux (Q.,) at -OWSeF, 
OWS N, and OWS V. The vertical dashed lines represent 
t+he mean of the distribution. 





Hebe 5-1, Statistical characteristics of the histograms of wind (u,), 
turbulent kinetic energy ca , and turbulent heat flux 
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Tables 5-2 to 5-4 present a summation of the monthly means, and 
standard deviations about these means, of the principal air-sea para- 
meters used to determine the atmospheric forcing depicted in the histo- 
grams. A smaller percentage of vapor pressure differences CE WE.) 
were calculated (see note 2) because atmospheric moisture information 
was not always recorded. 

The marine winds Cu.) are strongest at OWS P and weakest at OWS N, 
which accounts for the relative magnitudes of the uy, and ne distribu- 
tions. However a more significant observation is that at OWS P and 
OWS V the increase in The magnitude of the winds, from September To 
December, is nearly double that observed at OWS N. The air-sea tempera- 
Ture difference oe ae) and saturated vapor pressure difference CE -E.) 
are the additional parameters that are important for determining the 
turbulent heat fluxes. Tables 5-2 to 5-4 show that both of these para- 
meters are significantly larger at OWS V and OWS N than at OWS P. As 
@ consequence the smallest turbulent heat fluxes occur at OWS P, despite 
the largest observable winds. 

The magnitudes of these air-sea differences are related to the loca- 
Tion of the three stations relative to the mean atmospheric circulation 
in the North Pacific. The general westerly flow at OWS V and northeasterly 
Tlow at OWS N, together with the subsidence associated with the subtropical 
high pressure belt, continuously brings cold, relatively dry air in con- 
tact with the warm ocean at these two stations. At OWS P, however, the 
mean flow is westerly, the air mass has had considerable contact with the 
underlying ocean, and consequently the air-sea temperature and vapor pres- 
Sure differences are smal!. Additionally it should be observed that these 
differences are much more variable at OWS P and V than at OWS N (compare 


the means and standard deviations at each station in Table 5-1). This 


eS 


TABLE 5-2. Variability of atmospheric and oceanic parameters at 


OWS P. 


VARIABLE SEPT 

u_ (m/sec) oa 
= (4.6) 
T -T (°c) ce 
eS Gles0)) 
Eee (mb) i 
oe (1.8) 


Ore: 


OCT 


Ane) 
C565) 


NOV 


168 
CS) 


DEC 


eo 
CGe6) 


(1) First value represents the monthly/seasonal mean of the 
three-hourly readings while the second (in parenthesis) 


is the standard deviation. 


(2) Statistics are based upon 90% (21,030) of the possible 
three-hourly observations of E -E_ and 96% (25,524) of 
the remainder of the variables” “® 





TABLE 5-3. Variability of atmospheric and oceanic parameters at 


ows V. 

VARTABLE SEPT’ OCT NOV DEC 

iecin/ sec ) 6.7 7.8 G2 arora 
4 (3.7) (3.8) @aeth) (4.8) 
ieet (°C) eo le 220 Dei 
eo (1.4) (1.6) (2.0) (29) 
E -E (mb) 8.1 We 8.7 8.0 
~ 3 (i2) (4.4) (4.3) Cue) 


NOTE: 
(1) First value represents the monthly/seasonal mean of the 
three-hourly readings while the second (in parenthesis) 
1s the standard deviation. 


(2) Statistics are based upon 79% (10,930) of the possible 


three-hourly observations of E -E_ and 96% (14,065) of 
the remainder of the variables. 
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TABLE 5-4. Variability of atmospheric and oceanic parumeters at 
OWS N. 
VARTABLE SEPT cer NOV DEC 
u_ (m/sec) 5.6 5.9 6.9 7.3 
? (a9) 2.9) (3.3) (3.4) 
mee (°C) 0.9 1.2 7S res 
i (1.0) IRD) (1.3) Gia) 
E -E (mb) ee 8.4 8.2 i? 
4 Oe) O27 (2.9) (3-0) 
Pork: 


el) 


(2) 


First value represents the tmnonthly/seasonal mean of the 
three-hourly readings while the second (in parenthesis) 
1s the standard deviation. 


Statistics are based upon 76% (17,146) of the possible 


three-hourly observations of E “EL ands dS oe (25.57 7 aor 
the remainder of the variables. 
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would explain the additional variance in es at these two 
stations. 

The magnitude and variability in these data suggest strongly that 
the characteristics of the marine atmosphere at OWS P and OWS V are 
closely related to the frequency and intensity of the large winter 
stoms that occur at these locations. The atmospheric forcing at OWS N, 
on the other hand, is more closely coupled to the properties of the mean 
eiroularion. 

To examine the relative distribution of the turbulent kinetic energy 
and the turbulent heat fluxes, the percentage of these quantities along 
with the percentage of observations in each class interval were computed. 


For example 


nN 
3 
PE, = 100 Dt > ee (5-1) 


represents the percentage of turbulent kinetic energy That occurs in the 
a interval, n is the number of observations in the interval, and N is 
the total number of observations. Percentages of The turbulent heat 
fluxes (Q.) and wind (u,) were calculated with similar expressions. These 
percentages were accumulated from the smallest to the largest values and 
are compared in Fig. 5-3 in the form of cumulative frequency diagrams. 

The values for the cumulative percentage of lige and Q. are determined from 
The upper scale of the abscissa, while the cumulative percentage of obser- 
vations are obtained from the lower scale. For example, Fig. 5-3A indi- 


cates that all observations of 1 (ess Shaan 10° ony sees account for 


74%/90%/97% of the total observations at OWS P/V/N, but account for only 
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Figure 5-3. Cumulative percentage of turbulent kinetic energy (TKE) 
surface cooling, and observation at OWS P, OWS N, and 
OWS V. (See text for explanation.) 


I 


118 


‘= @ = — ; 
im Vi” ) ae 
| - 
ah 
Ss a <a eien’ 
OM aE ar ig 
SoA as 
(7712+) 4 
ore or i 





at fi at 
bai “at's j j 


a a 9 SF YY 
P cae 2 ae ee 
/ 
| 





30%/56%/83% of the turbulent kinetic energy flux computed from these re- 
cords. Similar comparisons can be made regarding the distribution of 

the turbulent heat fluxes from Fig. 5-3B (except that increasing values 
are read to the left). The horizontal displacement in these curves indi- 
Gates that the turbulent kinetic energy exchange at the three stations 
increases from a minimum at OWS N to a maximum at OWS P, while the turbu- 
lent heat fluxes are smallest at OWS P and largest at OWS V. The example 
just presented indicates that at all-stations a large percentage of the 
observations account for a relatively smaller percentage of the turbulent 
Kinetic energy fluxes. 

It is easily verified from Fig. 5=-2B that this is also valid for the 
turbulent heat fluxes. Furthermore the characteristic shape of the u,° 
curves suggest that although the magnitude of the turbulent kinetic energy 
flux is different at each stationits distribution is quite similar. The 
Shape of the os ae would indicate that the distributions of turbulent 
heat fluxes are similar at OWS V and OWS P but slightly different at 
OWS N. 

Additional information may be obtained by plotting the cumulative per- 
centages of U's a and uy, against the cumulative percentages of observa-~ 
Tions, as in Fig. 5-4. The first important observation from these curves 
is the relative invariance in the diStribution of u, and ae at the 
three stations (in these figures there is less than 3% scatter). Although 
These curves were constructed from the entire sample of values, curves 
drawn from the individual monthly values were nearly identical, and indi- 
cate That these distributions are also invariant with respect to month 
during the cooling season. The distribution of 0. was Identical to u, at 


OWS N while at OWS P and OWS V the curves indicate a higher percentage of 
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Figure 5-4. Percentage of ee oa as a function of percentage 
of observations. 
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large turbulent heat fluxes. For example 60% of the observations ac- 
count for 40% of the cooling at OWS N but only 35% at OWS P and OWS V. 
Again this is due to the additional variance tntroduced by the air-sea 
temperature and vapor pressure differences at OWS P and OWS V. 

An interesting observation from the ie curve is that it is repre- 
sented by the function ae = 200 PE - PEA » where FE is the cumulative 
percentage of observations and PE is the cumulative percentage of Te 
This function represents the locus of the circle centered at (0,100) with 
a radius equal to 100. A simple geophysical interpretation of this pro- 
perty ts that a considerable number of smal! values of the wind account 
for only a smal! percentage of the total turbulent kinetic energy flux, 
while a few large values represent a considerable percentage of the total 
flux. Additionally, but to a lesser extreme, the oe curves demonstrate 
a similar relationship for the turbulent heat fluxes. 

These relationships suggest strongly that a significant percentage 
of the turbulent kinetic energy and surface heat flux is exchanged at the 
air-sea interface during relatively short, but strong, atmospheric forc- 
ing events. Furthermore it is postulated that these events are directly 
assoctated with the eee of extratropical cyclones at OWS P and OWS V. 
At OWS N, however, these events are probably related to a pulsing in 
The mean flow, resulting from an alternate strengthening and weakening 
of the north-south pressure gradient as these storms pass to the north 
of OWS N. Additional evidence to support this assumption will be pre- 
sented in the final section of this chapter when the distinguishing 


emeracteristics of these events are examined. 
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C. COMPARISON OF THE LONG-TERM MEAN FORCING AND 
BEEANIC RESPONSE AT THE OCEAN WEATHER STATIONS 


In this section the long-term mean daily averages, depicted previous- 
ly in Fig. 3-l, will be examined to determine the mean energetics of the 
forcing at the three ocean weather stations and the magnitude of the 
mean oceanic thermal response. The quantities will be used in subsequent 
analysis as a basis for Josoamining what constitutes strong forcing and 
strong oceanic response. We will also present additional evidence from 
these data that is supportive of the assumption that the character of 
the atmospheric forcing at OWS P and OWS V is determined by the synoptic 
Storm patterns, and at OWS N by the mean flow. Finally simple one-dimen- 
Sional reasoning will be applied to these data to demonstrate that the 
major features of the long-term mean thermal response are explainable in 
terms of simple vertical processes. 

To obtain the energetics of the forcing, the daily mean values of 
i, Qo? Qe and OF were summed and the monthly totals are presented in 
The first four columns of Tables 5-5 to 5-7. The next two columns show 
the long-term mean monthly changes in the sea-surface temperature (ASST) 
and the mixed layer depth (AMLD). The relative contribution of the net 
surface heat flux (20.) to the mean sea-surface temperature change is 
presented in the last column as ASST(2Q_). This contribution was esti- 
mated by distributing 10, over the mean layer depth at the beginning of 
The month (upper value), and at the end of the month (lower value). 

Since the layer deepens between these two depths, the actual contribution 
of 20, is somewhere between these two extremes. 

The seasonal trends in the forcing may be examined to show that both 
The magnitude and variability of these quantities correlates well with 


The locations of the ocean weather ships relative to the major storm 


eZ 











9° 0= (L°0) (9TE°O) (9€0°0) (hHEO°O) (EET °O) 
L°0- LT- 2°T- ero O= 9T0°0 LL0°0- GLE°O ext! 
(6°0) (L8T°O) (9S0°0) (e200) (6hT°0) ae! 
tz- | O°2- H9S°0- e200 T80°0- L8€°0 \ON 
Cet) (O6T°O) (TT0°0O) ($Z0°0) Gror0) 
he- 8°Z- GZE°0- LhO°O 640°0- See LOO 
(680°0) (4T0°0O) (Gn0G0) (Z90°0) 
OT- 69T°0 €90°0 oT0n 0 ELE 0 Idd$ 
(OW/D,) | (OW/W)} (OW/D,)}) CULUoW/AT) | = CyZUoW/AT) | (YZUOW/AT) (6098/ UO) | HINON 
(“OZ)LSSV CIWV 
(*uoTeUeTdXS JOJ 7x92 
28S) ‘"ueew ATYUOU ey ANOGe UOCTSTASeP PXeEpueAS 3} ST STSayjUeTed UT SNTeA 
*(adertene JeeK 47) g SMO 3e esucdsei TJeAePT pexTu pue BSuUTOIJOJ uPoUul uIa}-sUuOT =°G-G TIAVL 


on OF 
Us 






u 5S e 3 
C01) Oz (,_0TX) OF (,_0T%) Oz (57 -0T*) ny 


123 





olla cD COnHE"O) (S00 °0) CZe0) 0) (THO*O) - 
oc s Cle SC nee ie 8h0°O Sena CTé°0 JAG 


Tas (8°T) (962 °0) (400°0) (TE0°O) (9hO°O) 
ae EGe GG 906°0- 8S0°0 641° O- 09T'0 AON 





8°0- Chae (S8Z°0) (800°0) (9Z0°0) (8€0°0) 

Z°T- LT- | T'?- 90h*0- L80°0 LZU°0- HOT’ O 100 
Cae CEB) (812 °0) (600°0) Cocos (ZHO°O) ; 
€°O+ 0°T- 0°T- SL0°0 Hee 90T°0- TOO ids 


(OW/D—) | (OW/W) | (OW/D5)} CufUoW/A]) | (YfuoW/A]) (YfuoU/A | ) 


Cozyussv | any | assy (, 01x) 0g (,_0TX) "03 (._0TX) "02 ( 


Bean HLNOW 


ENE 
2-0 de 4 


(*uoT eueTdxXe JOJ 21x93 
38S) ‘ueow ATYZUOW BY ANOGe UOTSTASP PrepueS sy} ST STSeyZUeTed UT SNTeEA 


‘(adeuaae week sT) A SMO 2e esuodseu UehePT pextw pue BuTOUOy ueaW uIE}-BUC] «*9-S TIAVL 


12:4 








Uu 
( 03)LSSV CIWV 1Soe 


IH0= (2°T) (9h? °0) (S00°0) ($20°0) (£20°0) 
CUS ule La bea Us 6S0°0 ELLOS 8Z0°0 JAd 


9°0- Caz0)) (SEZ°0) (400°0) (eculo) (6Z0°0) , 
8°O- le C°1- 08h’ O- 0G STT°0- S90°0 AOQh 





6 Us (6°0) CESie0) (Z00°0) Con <0) (ST0°O) 
ces ols 6°0- ies 960 °0 80T°O- THO’ O TOO 






9°0+ (6°0) (96T°0) (0T0°0O) (HT0°0O) (9€0°0) 
8°O+ TI- 0 nt C0 6LL'0 960°0- SE0°0 12s 
(OW/D9) | (OW/W) , (CW/O,} (UsUOW/A] ) (Y4UoUW/A | ) (yfuow/Ay ) Opes et) HLNOW 


ue a= tok : 3 
Cor) O% au) OZ CEs) Oz (7,01), nz 


(*uoT}ePUeTUXe UOF 1X97 
39S) ‘wesw Apy.UOW BY} JNoqe UOTSTASP PAepueZS sy} ST STSeYyZUSeTed UL SsNnTeA 
‘(aseraae Teed ¢7) N SMO 3e esucdsad sweAPT PaxTW pue BuTOUOJ uPA WIe-BSUOT °*1-S TIYVL 


mn eet eee” ene a 


2 


> 
} 
- > 
= 
_- 
a 
—_ 
- 
- 
> 
; ws 
. 
- 
- 
oe ee 
= 
- 
— 





tracks in the North Pacific. As expected, the largest flux of turbulent 
kinetic energy Oe is received at OWS P and the smallest at OWS N, 
with a general increase in magnitude during the season at all stations. 
At OWS P and OWS V, however, ri Increases rather abruptly in October 
and November respectively. It is quite evident that at OWS N no abrupt 
changes in Ble occur during the cooling-season in any month. This is 
rather suggestive that at OWS P and OWS V these abrupt increases in ee 
mark the onset of the winter storm season at these stations. At OWS N 
the trend in re would suggest a slow intensification of the wind, which 
would be indicative of the wintertime strengthening of the mean circula- 
ion. 

The monthly Q, values at the three stations show the same relation- 
ships as previously noted in the histograms (Fig. 5-2). The most intense 
Turbulent heat fluxes occur at OWS V and the weakest at OWS P. By compar- 
ing the mean and standard deviations of 20. at the three stations, it may 
be noted again that the turbulent heat fluxes are more variable at OWS P 
and OWS V than at OWS N. As expected, the effective solar insolation, 
LQ.» decreases throughout the season and with latitude. Additionally it 
Should be observed that the insolation at OWS P has much more variability 
Than at the other stations. Since this can only be introduced by the 
total cloud cover (see Appendix A), it is indicative of the changes in 
the cloud patterns that accompany the cyclonic storms at OWS P. At OWS P 
The cloud-types are normally low stratus and the sky frequently covered 
by heavy overcast throughout this season. The exception to this pattern 
Occurs in the cold, dry air mass behind the cold front where the cloud 
patterns are broken and, therefore, allow more insolation to reach The 


sea surface. At OWS V and OWS N, however, the sky is not normally overcast 
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and the storm systems do not introduce as much variability in the total 
cloud amounts. Finally the net heat gain or loss by the ocean is re- 
flected in the fourth column of Tables 5-5 to 5-7 as xQ0° There is a 
net heat gain, at all stations, during September and this downward heat 
flux is largest at OWS N and smallest at OWS V. During the remaining 
Three months 20, is increasingly negative, at all stations, and is most 
negative at OWS V and least negative at OWS N. 

The oceanic thermal response has characteristic trends, at each sta- 
Tion, that we will attempt to explain with one-dimensional reasoning. 

The total mixed layer depth change (AMLD) is largest at OWS P and smallest 
at OWS N; however, the trend is similar at all three stations. The deep- 
ening rate increases early in the season and is reduced at the end. The 
seasonal sea~-surface Late change (ASST) is comparable at OWS V and 
OWS P, but significantly smaller at OWS N. It should be noted that the 
trend in ASST is different at OWS P than at OWS V and OWS N. At OWS V 

and OWS N, ASST increases throughout the season, while at OWS P ASST in- 
creases during October and decreases in the final two months. 

If one-dimensional reasoning is applied to the long-term forcing and 
oceanic response, a considerable amount of the observed changes are ex- 
Platnable. Consider first the amount of the monthly sea-surface tempera- 
Ture change (ASST) that may be explained by the net surface Tlux (LQ). 
september, at all stations, marks the transition between the heating and 
cooling seasons and the data indicate that the surface fluxes increase 
the heat content of the upper ocean at all stations. This occurs during 
a period when the sea-surface temperature is normally decreasing at OWS P 
and OWS V, and unchanging at OWS N. Therefore the sea-surface tempera- 
Ture change, during September, cannot be explained (and certainly not 


Predicted) in terms of the surface heat fluxes alone. The net surface 
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fluxes become increasingly negative from October to December; however, 
ASST(2Q_) increases from October to November and then decreases in Decem- 
ber. This is because the ratio of the net surface fluxes to the mixed 
layer depth increases from October to November and decreases in December 
ie. ASST (20) = = 9). This illustrates the importance of the mixed 
layer depth in determining the effect of surface fluxes on the sea- 
surface temperature HOF os: 

The relative importance of the surface fluxes to the total sea-sur- 
face temperature change may be ascertained by comparing the cumulative 
Bemttibution of ASST (20? as Wie Cet S22 Ge B75 ele siugies 1 ime, eats rete: 
months, with the total sea-surface temperature change during these months. 
At OWS V and OWS N greater than one-half of ASST may be directly explain- 
able by ASST (LQ), while at OWS P only one-third is explainable. Final- 
ly, aS a measure of the dominance of the surface fluxes, consider the 
months during which ASST (2Q_) accounts for greater than 50% of ASST. At 
OWS V and OWS N these fluxes become dominant during November while at OWS 
P this does not occur until December. It is evident that, even during 
The months when these surface fluxes appear dominant, they are not capable 
of explaining the total sea-surface temperature change. 

The other significant vertical process which might explain the differ- 
ence between ASST Bae Soi (2Q,) is the heat flux at the base of the mixed 
layer that occurs during entrainment. In fact, during large deepening 
events, the entrainment heat flux may actually dominate the surface heat 
flux, as was demonstrated in the previous chapter. Since this heat flux 
Cannot be computed from these data,we must infer its relative magni tude 
from the turbulent kinetic energy flux eae the mixed layer depth 
Change (AMLD), and the characteristic thermal structure at the three ocean 


weather stations. 
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It is important to observe that nearly twice the turbulent kinetic 
energy is transferred to the ocean at OWS P as at OWS V, and greater 
Than five times as much as at OWS N. Additionally, the maximum seasonal] 
residuals between ASST and ASST (20? that must be explained are -5.2°C, 
-3.7°C, and -2.4°C at OWS P, OWS V, and OWS N respectively. The rela- 
tive magnitudes of these fluxes and residuals are understandable in 
terms of (2-49), (2-50), and the fundamental! differences in the thermal! 
structure at the three stations, illustrated in Fig. 5-5. The character- 
istic shape of the profiles at OWS V and OWS N are quite similar, but 
considerably different than the profile at OWS P. This is because at 
both OWS V and OWS N the parent water mass is North Pacific Central 
(Sverdrup, 1942), while the water mass at OWS P is Pacific Subarctic. 
However, the important difference to be observed in the three profiles 
is the magnitude of the temperature gradient immediately below the mixed 
layer. The temperature change from /0-80 meters is approximately 2.5°C 
/2.1°C/1.5°C at OWS P/OWS V/OWS N. The relative magnitudes of the 
temperature difference (AT) at the base of the mixed layer is in The 
same proportion at the beginning of September. 

Equations (2-49) and (2-50) indicate that the largest amount of tur- 
bulent kinetic energy would be required to mix the upper layers of the 
ocean at OWS P and the least at OWS N. This is completely consistent 
with the relative magnitudes of Bue noted previously in Tables 5-5 To 
2-7. Additionally, comparing AMLD and AT at each station, we should ex- 
pect the largest entrainment heat flux (SP AT) at OWS P and the smallest 
at OWS N. This is consistent with the relative magnitudes of the sea- 
surface temperature residuals computed at the three stations. 


The seasonal trend in the mixed-layer depth is understandable in 


terms of the magnitude of ee and the concept of dissipation enhancement, 
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Figure 5-5. Typical mid-season temperature profiles at the three 
ocean weather stations. 
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discussed in the previous chapter. As The layer deepens, the entrainment 
zone is displaced from the surface production zone, and a greater percent- 
age of the turbulent kinetic energy is dissipated. In accordance with 
(2-33) as the ratio of G,-D, and h decreases the mixed layer deepening 
rate also decreases. This occurs in November at OWS P and OWS N and in 
December at OWS V. 

lt is encouraging to note the major trends in these data are consis- 
tent with simple one-dimensional reasoning. However, it is important to 
realize that these trends are not explainable solely in terms of the 
Surface forcing. It is necessary to be able to specify both the surface 
heat flux and the entrainment heat flux to accurately describe the sea- 
sonal mixed layer evolution. The entrainment heat flux will dominate 
The surface heat flux early in the season when = and AT are large and 
The surface heat fluxes are small. By the end of the season the large 
Surface heat fluxes become the dominant factor in changing the sea- 
Surface temperature. 
See MEORTANCE OF LARGE ATMOSPHERIC FORCING EVENTS IN 

DETERMINING THE SEASONAL EVOLUTION OF THE MIXED LAYER 

In this section the seasonal evolution of the mixed layer will be ex- 
amined relative to the long-term mean trend to illustrate the importance 
of large atmospheric forcing events in changing the upper ocean thermal 
structure. The purpose of this analysis will be to demonstrate that the 
Timing of the first large storms (whether they occur early or late in the 
season) is very important to the seasonal evolution of the mixed layer. 
Evidence will be presented to illustrate that the underlying thermal 
structure is as important as the surface forcing. Two cooling seasons, 


Taken from the OWS P record, were chosen to demonstrate these principles, 
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because both the atmospheric forcing and the thermal structure were 
very different during these two years. 

Mid-October temperature profiles (Fig. 5-6) illustrate that the mixed 
layer was deeper and cooler during 1963 than during 1959. Notice also 
that the temperature gradient at the base of the layer during 19635 was 
considerably weaker than during 1959. The final observation is that in 
1963 there was a much stronger thermal gradient below 60 meters during 
Pee than in 1959. 

Figures 5-7 and 5-8 depict the observed forcing and oceanic response 
during these two seasons, with the long-term mean trends superimposed. 

In Table 5-8 the strength of the forcing and response characteristics 
are presented and compared with the long-term mean. In Figs. 5-7A,B it 
is observed that at the start of the 1959 season the mixed layer was 
considerably deeper and cooler than the long-term mean. However, at 
The beginning of September 1963 a very warm, shallow layer is present 
eeigs. 2-8A,B). 

The character and timing of the atmospheric forcing during these two 
seasons may be observed in Figs. 5-/7C,D and 5-8C,D. During 1965, the 
period from 15 October to | November marks the beginning of the winter 
storm season. In 1959, with the exception of the one significant event 
in mid-October, the large forcing events do not begin until mid-November. 
The difference in these storm patterns coupled with the different initial 
temperature structures result in the evolution of the upper ocean being 
quite different during these two seasons (compare Figs. 5-/7A,B and 5-8A,B). 
In the previous chapter, it was demonstrated that large forcing events 
that occur early in the season produce much larger oceanic thermal re- 
sponses than events occurring late in the season. This is consistent 


with what is observed in Figs. 5-7 and 5-8, and we will examine the 
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Atmospheric forcing and oceanic response at OWS P for 


the 1959 cooling season. 
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Figure 5-8. Same as Figure 5-7 except for the 1963 cooling season. 
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TABLE 5-8. Atmospheric forcing and response at OWS P for the years 
aoe. o6s- 
Forcing and Response (as a percentage of the 
seasonal total) 
1963 1959 
Month ‘%u,, oQ, Q. 2Q, %AMLT SAMLD Month %u,, Qa, Qa QQ, SAMLT %SAMLD 
Sept 12 20 42 -11 25 21 Sept ik 14 42-12 -6 -16 
Det 40 4O 34 4 ili SOGIE OA 2 ee SO Ae aq Ze. 
Nov yee 26 15 Hi 16 “ Nov Sop 2516 33 Mise 
Dec Zo 14 Seco 8 Dec Zo 0) a a oe 41 
“Missing Data 
Foreing and Response (as a percentage of the 
long-term mean) 
1963 1959 
Month %u,,” %Q, %Q. %Q, SAMLT %AMLD Month %u,,” 3Q, %Q_ %Q, %AML %AMLD 
Sept 95 134119 8. 248 150 Sept 115 92 99 11€ -28  -80 
Oct imo 156 130 1- dey LOG = Oet 83 seme 35" bo ot 98 
Nov ee 96 10/7 EB = Nov 121 iat ©6=— 8 
Dec 32 Doster 4. AE “ bee sot Se aero wisi. abe) 
season 1i4 108 120 9h 37 100 Season 1C& 09 99 121 #765 dal 
Oceanic Response (.36?) Oceanic Response (1959) 
Month — MLT(“C) MLD(M) Month MLT(°C) MLD(M) 
Sept 5.3-12.9 be 33 sept Zoe .6 45-37 
Set 2. 9- 8.1 -/9 Wer 12.8-11.4 37-51 
Nov Se 6.5 79 Nov 11.4- 9.4 51-* 
Dec 6.6- 5.8 -30 Dee C.4- 7.2 * 96 
Oceanic Response 
(long-term mean) 
Month — MLT(°C) MLD(M) 
2 PSH 1263 28-38 
aie P2365 38-62 
Nov Seo= 5 62-83 
DEL Leo 6. 3 83-100 
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oceanic evolution during these two seasons To show that the major features 
are explainable by one-dimensional reasoning. 

During 1963 the month of September had weaker than normal downward 
fluxes of turbulent kinetic energy (uy?) and heat (Q.) (compare percentage 
of long-term mean in Table 5-8). The mixed layer, however, deepened and 
cooled more than the long-term mean. The larger than normal deepening and 
cooling rate was the result of the layer being warm and very shallow at 
the beginning of the period. Because the layer was shallow, only a rela- 
tively small amount of turbulent kinetic energy was dissipated, which 
resulted in the large entrainment rate. As would be expected, the verti- 
cal re-distribution of the heat in this warm layer resulted in a larger 
Than normal entrainment heat flux and consequently abnormal sea-surface 
Sererature change. During October 1963 the large winter storms are 
responsible for an energy exchange between the atmosphere and ocean which 
is significantly larger than the long-term mean. These large fluxes force 
The mixed layer to deepen and cool, such that by the end of the month the 
layer is deeper (8 m) and cooler (1.4°C) than the long-term mean. During 
these first two months 52%(60%) of the seasonal ug (Q,) was received by 
the ocean, which is 10% greater than the long-term average. However, this 
abnormal forcing accounted for 75% of the mixed layer response during 
the season, which is 25% greater than normal. This large response occurred 
because the strong forcing came early in a season when a large heat 
storage had taken place in a very shallow mixed layer. 

The final two months of 1963 were characterized by weaker than normal 
forcing. The upward heat fluxes (Q,) account for only 40% of the seasonal 
total (10% less than normal) and this is a consequence of the anomalously 
cool sea-surface temperatures. The turbulent kinetic energy flux (uy?) 


received during the final two months jis less than the long-term average 
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but is actually about the same amount received during the first two 
months (48% of the seasonal total). However, the response of the layer 
is only 25% of the seasonal total. The turbulent kinetic energy flux 
received during the final two months is simply not sufficient to accom- 
plish any significant deepening. Evidently this ts because the entrain- 
ment zone is deeper than normal, a larger percentage of the Turbulent 
kinetic energy is dissipated, and the ratio of G,-D, to h is less than 
normal. 

Examination of the total seasonal forcing (Table 5-8) shows that 
greater turbulent kinetic energy and smaller net heat fluxes were ex- 
changed at the air-sea interface during the 1963 season, which resulted 
in a normal seasonal deepening rate but a larger-than-norma| sea-surface 
temperature change. This is consistent with our earlier findings that 
a temperature structure that has a large heat storage in the near surface 
layers will require a large amount of energy flux to deepen the layer. 
Additionally this deepening will be accompanied by a large entrainment 
heat flux and, consequently, a large sea-surface temperature change. 

This example demonstrates clearly that the energy fluxes received 
early in the season results in a much larger mixed layer response than 
comparable energy fluxes received late in the season. An additional ob- 
servation is that during September and October the anomalously large up- 
ward heat fluxes (Q,) are accompanied by significantly larger-than-norma | 
downward fluxes of solar radiation (Q.)- This supports our previous 
Suggestion that the large variance in Q. (Table 9-5) at OWS P is due to 
the large winter storms. The air mass is generally cold and dry during 
stormy periods and there is usually less cloud cover than the mean winter 
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The evolution of the mixed layer during the 1959 season was quite 
different than during 1963 (Figs. 5-8A,B). Since the mixed layer was 
deeper and cooler Than normal at the start of the season, if is pre- 
Sumed that a significant re-distribution of the heat was accomp! ished 
by strong summer forcing (evidence will be presented to support This 
presumption). The Turbulent kinetic energy was stronger than normal 
during September 1959, but so was the downward net surface heat flux 
(see Table 5-8). As a result, the. mechanical energy was insufficient 
to mix this additional heat to the old mixed layer depth and a new 
shallower mixed layer formed above the old one (see Figs. 5-8A,B). By 
the end of the month the layer depth had retreated to near The long- 
term mean while the temperature is warmer than normal. The forcing re- 
ceived during October 1959 is much weaker than normal and only a weak 
mixed layer response takes place. The forcing received during the month 
of November is characterized by extremely strong mechanical forcing 
while the upward heat fluxes are smaller Than normal. Because of the 
gap in the bathythermograph record it was not possible to determine the 
mixed layer depth at the end of The month. However, Fig. 5-8B would 
seem to indicate that the deepening rate is larger than normal. Because 
The net surface flux (05) is less Than normal the sea-surface tempera- 
Ture decreases only a normal amount and at the end of the month the 
Aayer is probably deeper but significantly warmer than normal. December 
1959 is characterized by average mechanical forcing but extremely larger 
upward turbulent heat fluxes. The sea-surface Temperature change is 
much larger than the long-term mean as a result of the large surface 
heat fluxes. At the end of the season the layer depth is about normal 


but the sea-surface temperature is approximately |°C warmer than normal. 
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The seasonal totals in Table 5-8 show that the forcing was actually 
stronger than the long-term mean and yet less than normal oceanic re- 
sponse occurred. This is primarily the result of the deep mixed layer 
that was present at the beginning of the season. The forcing during the 
first two months was not sufficient to establish the normal deepening 
and cooling rate. The large winter storms came late in the season and 
were not able to reduce the sea-surface temperature the additional amount 
necessary to bring the ocean back to The long-term mean. 

These two examples illustrate that the timing of the strong atmos- 
pheric forcing events coupled with the characteristics of the underlying 
thermal structure are of primary importance in determining the evolution 
of the upper ocean during the cooling season. Once again it should be 
noted that the large forcing events produce a much greater oceanic re- 
sponse when the mixed layer depth is shallow (usually early in the season). 

lf may be seen in Fig. 5-5 that the temperature structure is quite 
different at OWS P than at OWS N and OWS V. Furthermore Fig. 5-6 shows 
That the thermal structure may have considerable variability from year to 
year at a particular location. Since the seasonal thermocline is estab- 
lished during the spring and summer, these figures would suggest that 
the mechanisms responsible for its formation may vary depending on loca- 
tion and year. The spring and summer forcing, observed during 1959 and 
1964 at OWS P, were compared to see if the differences in the temperature 
structure, depicted in Figs. 5-6, 5-7A,B, and 5-8A,B, are explainable 
in terms of one-dimensional processes. The mechanical forcing observed 
during the spring in 1959 was much weaker than observed in 1963. However, 
the synoptic patterns observed during the summer indicates that the 1959 


season was characterized by a high incidence of strong storm activity 
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while the forcing was extremely weak during 1963. Consequently the 
layer retreated very rapidly during the spring of 1959 and because very 
little heat was mixed below 60 meters, the seasonal thermocline below 
this level was nearly isothermal. On the other hand, during 1963 a 
considerable amount of heat was mixed into the deep layers by the strong 
spring forcing, and a relatively strong thermal gradient was established 
below 60 meters (see Fig. 5-6). The strong summer forcing during 1959 
caused a large downward transfer epeneat and resulted in a very deep, 
cool mixed layer at the beginning of September (see Fig. 5-7A,B). The 
weak forcing during the summer of 1963 was insufficient to maintain a 
normal summertime mixed layer depth and a very warm, shallow layer 
developed (see Fig. 5-8A,B). 

This simple analysis suggests that the synoptic storms may be as 
important to the establishment of the seasonal thermocline during the 
heating season as they are to the subsequent erosion that takes place 
during the fall and winter cooling season. The evclution of the upper 
ocean during the cooling season is determined to a large extent by the 
nature of the stability of the thermocline. Therefore it is important 
to realize that the entrainment process may be influenced by the nature 
of the forcing during the previous spring. A simple conclusion would 
be that the spring and summer heating season should be examined in 
detail to determine the relationship between the principal mechanisms 
that govern the formation of the seasonal thermocline. 

E. IMPORTANCE OF THE STRONG ATMOSPHERIC FORCING 
EVENTS AT THE THREE OCEAN WEATHER STATIONS 
In this final section the total record (see Table 3-|) of the surface 


and near-surface parameters will! be examined to determine the relative 


141 


r 





importance of large forcing events to the total evolution of the upper 
ocean thermal structure. The objectives of this analysis will be To 
determine: 


(1) the principal characteristics of the large events at the Three 
Siawions; 


(2) the percentage of the total energy that is exchanged during 
These events. 


Figure 5-9 is presented to explain how The analysis was performed. 
An event will be defined as those individual periods when the forcing 
was greater than the long-term mean. For example, in Fig. 5-9 there are 
12 individual events during the season. If thé forcing is greater than 
The long-term mean at the beginning or end of the season it is counted 
as an event. If an event takes place during two months (for example 
events 3 and 9), the amount of energy received during each month is 
calculated and the event is credited to the month in which the largest 
energy exchange occurs. For example the third event in Fig. 5-9 would 
be credited to October and the ninth event to November. The duration 
of an event is defined as the time during which the forcing is greater 
Than the long-term mean. An additional parameter was the ratio of The 
peak value of the forcing to The long-term daily mean. This peak-to- 
mean ratio (Pk) was used to examine the effects of the larger events. 
Soeeinstance we might only examine events where Pk > |.5, etc. The 
analysis was performed first by defining the forcing events in Terms of 
ne m7 curves (mechanical events). Next the analysis was repeated 
with the events defined by the Q. curves (cooling events). The results 
from both of these definitions will be presented. 

The response of the ocean thermal structure was estimated from the 


seasonal sea-surface temperature change (based on the surface marine 
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Figure 5~9, Characteristic of events. 
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observations) that occurred during each event. The ratio of the cumula- 
tive changes in sea-surface temperature taking place during the events 

to the total seasonal change is the percentage of the total response 

that occurs during the events. In a similar manner, the percentage of 
the total duration, ee and om occurring during these events was calcu- 
lated. The large data gaps in the bathythermograph record would not per- 
mit analysis of mixed layer depth during these events. However, the 
examples presented earlier in this chapter, and the results from the 
numerical models indicate that large changes in sea-surface temperature 
are normally accompanied by large changes in the mixed layer depth. 

In the first part of this analysis we will examine the percentage 
of the monthly forcing and response that takes place when the forcing 
Ceither nee or oe is greater than the long-term mean. The results of 
this analysis is presented in Tables 5-9 to 5-I|! and is based upon the 
24/23/15 seasonal records considered at OWS P/OWS N/OWS V. 

The data presented in these tables is consistent with our original 
hypothesis that a significant percentage of the atmospheric forcing and 
oceanic response takes place during periods when the forcing is strong. 
The first observation is the relative invariance in the duration of the 
forcing at all stations. The mechanical events characteristically are 
Shorter in duration than the cooling events because large air-sea tem- 
perature and vapor pressure differences persist longer than the strong 
winds. The relationship between the duration and the percentage of tur- 
bulent kinetic energy exchanged during these periods is also very consis-= 
tent at all stations. At all stations a large percentage of ee takes 
place during a relatively short time frame. Although the magnitudes of 


these large energy events are significantly different at the three 
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stations, these data would support the suggestion, made from the histo- 
grams, that the distribution of the is fairly invariant at the three 
Stations. Because OF includes both Q., and QO, it ts much more variable 
Than wee During September, for instance, when there is normally a net 
downward heat flux, a net upward heat flux takes place during these strong 
forcing periods (compare Oe at each station). This simply means That in 
September the surface cooling that does take place probably occurs during 
These strong forcing periods. It is also evident from these data That 
There is net heating taking place during October at all stations. This 
is reflected by 0, being greater than 100% during these months. AS ex- 
pected 7 is larger in the mechanical events and Q,, is larger in the 
cooling events. 

The final, and most significant observation, is that during all 
months (mechanical and cooling events) a significant percentage of the 
sea-surface temperature changes (ASST) takes place during a relatively 
short duration. Also the general trend in these data is that ASST de- 
creases as the seasons progress. This is true for both the mechanical 
and the cooling events. This once again supports the argument that the 
large atmospheric forcing produces a larger oceanic response early in 
The season when the mixed layer is shallow. 

The second part of the analysis was designed to determine the charac- 
teristics of these large events and their relative importance to the 
seasonal evolution of the thermal structure. The results of this 
analysis are presented in Tables 5-12 to 5-I4 in the order st imeres sinc 
peak-to-mean ratios (Pk). Therefore the statistics in Table 5-12 repre- 
sent every period in the record identified as an event, while those in 


Tables 5-13 and 5-14 represent only those events where Pe peed 
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Pk > 2.0 respectively. The first four entries in these tables (duration, 
re, Q0? ASST) are the percentage of the seasonal totals that occur dur- 
ing these events. The final three entries are useful for comparing the 
characteristics of the events at each station. 

The data in these tables supports the hypothesis that a significant 
percentage of the energy exchange takes place during these large atmos- 
pheric forcing events and results in a large oceanic thermal response. 
This may be easily verified by comparing the relationship between dura- 
maton, fe, Qn? and ASST for both the mechanical and cooling events listed 
in Table 5-12. Once again it may be observed from this table that, con- 
sidering all possible events, the percentage of duration of the cooling 
events is larger than the mechanical events. Comparing these same rela- 
tionships in Tables 5-13 and 5-14 shows the importance of The larger 
events to the upper ocean evolution. The general trend is for the larger 
mechanical events to remain significant while the larger cooling events 
become less important. For example the cooling events with Pk > 2.0 are 
almost insignificant. The exception to this trend is OWS V. IT should 
be observed in Tables 5-13 and 5-14 that for Pk > 1.5 ASST is only 31% 
pamens V, and for Pk > 2.0 SST is reduced to only 22%. These percentages 
are significantly smaller than at OWS P and OWS N. It was noted previous~ 
ly that the intensification of the atmospheric forcing does neon Ceeuls 
until November at OWS V. Therefore it is possible that the large events 
femee 1.5) occur later in the season at OWS V and this would account for 
the smaller sea-surface temperature response. Before discussing the 
Characteristics of the events one additional observation may be made by 
comparing uy? and Qn, during the mechanical events in the three tables. 


At OWS P and OWS V, the percentage of ie is consistently larger Than Qe 
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while at OWS N it is consistently smaller. This is additional evidence 
that the surface cooling at OWS N is more closely coupled to the wind at 
OWS N than at OWS P and OWS V. 

The final part of this analysis was to determine the principal charac- 
teristics of the large events. The final three quantities (number of 
events/month, duration of the events, and peak-to-mean ratio) will demon- 
strate that the character of the events is very similar at OWS P and 
OWS V but quite different at OWS N. In fact these quantities support the 
hypothesis that these events are directly related to the properties of 
the large winter storms at OWS P and OWS V, while they are related To 
changes in the mean circulation at OWS N. Table 5-I2 shows that the 
average number of mechanical and cooling events is largest at OWS P and 
smallest at OWS N. At all stations a larger number of mechanical events 
occur than cooling events but have a shorter duration and a larger peak- 
To-mean ratio. However it is important to observe that, for the mechani- 
cal events, the duration and Pk (mean and standard deviation) are very 
Similar at OWS P and OWS V and noticeably different at OWS N. This be- 
comes even more apparent for the larger events as may be seen in Tables 
9-13 and 5-14. The larger mechanical events at OWS P and OWS V appear to 
be well organized on time scales of 2-3 days, have similar Pk, and occur 
at the same frequency each month. These large events at OWS N however, 
eecur less frequently, are of longer duration, and exhibit a very large 
and variable peak-to-mean ratio. These are exactly the nature of the 
Statistics that would be expected if the events at OWS P and OWS V were 
directly associated with the passage of extratropical cyclone systems. 
The statistics at OWS N are compatible with the presumption that the 


events are related to a pulsing of the mean circulation. 
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The important implication from this analysis is that these large 
events occur frequently during the fall and winter at all stations. 
Furthermore they play a significant role in determining the character- 
istic evolution of the upper ocean thermal structure. Therefore they 
deserve special attention in any prediction scheme that attempts to 


reproduce this evolution. 
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VI. CONCLUSIONS 


The principal objective of this research was to examine the upper 
ocean thermal structure modifications that take place in response to 
strong atmospheric forcing events during the fall and early winter cool- 
ing seasons. The motivation was the fact that there was a serious gap 
in our fundamental knowledge regarding the role that these strong events 
play in the total thermal structure evolution during the fall and early 
winter, 

The one-dimensional hypothesis was evaluated at the three North 
Pacific Ocean weather stations using modified versions of the Kraus and 
emets (1967), Kim (1976), and Elsberry, ef al. (19/6) models. fhe 
Performance of the three models was evaluated using 49 independent data 
Sets from The historical series of marine observations at OWS P, OWS N, 
and OWS V. Additionally the EFT model was used to isolate and examine 
the relative importance of the principal one-dimensional mechanisms at 
the three stations. Finally, the large body of surface and near-surface 
marine observations were examined at the three North Pacific Ocean 
weather stations. A new analysis technique was employed to determine 
the relative importance of strong atmospheric forcing events in the Total 
fall and early winter thermal structure modifications at these stations. 
From the analysis the following significant conclusions were drawn. 

|. The integrated effect of the strong fall and winter atmospheric 
forcing events is the dominant factor in the modification of the upper 
ocean thermal structure at OWS P, OWS N, and OWS V. For example, these 


strong events occurred during approximately 35% of the time at the three 
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ocean weather stations. However, 85%/68%/57% of the sea-surface tempera- 
ture change at OWS N/OWS P/OWS V occurred during these periods. Observa- 
tions from the individual data sets indicate that one can expect similar 
responses for mixed layer depth changes. It was therefore concluded 

that these strong events can not be excluded from any forecast scheme 
developed to simulate upper ocean response during the cooling season. 

2. the response of the upper ocean during the strong events investi- 
gated in this study was largely one-dimensional. Additionally, a modi- 
fied version of the EFT model consistently demonstrated better agreement 
with observations than either the KT or KIM models. The KT model consis- 
tently predicted mixed-layer temperature and depth changes which were 
much larger than the observations. This result suggests that the fraction 
of turbulent kinetic energy that is available for entrainment is not a 
constant fraction of The turbulent kinetic energy transferred to the 
ocean by the wind, as postulated by Turner (1969). The superior perform- 
ance of the EFT and KIM models relative to the KT model suggests that the 
amount of wind-generated turbulent kinetic energy arriving at the mixed 
layer interface decreases as the mixed-layer depth increases. However, 
the exponential parameterization of dissipation enhancement employed in 
the EFT model was found to be more effective in preventing excessive 
deepening rates than the linear representation in the KIM model. IT was 
further concluded that the convectivel y-generated turbulent kinetic energy 
is largely non-penetrative, in agreement with Gill and Turner (1976). 

5. During these strong forcing events a large component of the sea- 
Surface temperature change is due to the vertical fluxes of heat at the 
surface and at the base of the deepening mixed layer. Depending on the 


magnitude of the turbulent fluxes of kinetic energy, the heat exchanged 


156 







DUMP Disc 


CLEAZ Sut dE 


¢ 
AcE 
oe a! =): hold 





at the air-sea interface, and the depth of the mixed layer, either of 
these heat fluxes may dominate the local heat budget of the mixed layer. 
However the important conclusion that may be drawn from this research is 
that an accurate specification of both fluxes is essential to understand- 
ing and predicting the sea-surface temperature changes. The effect of 
entrainment mixing is most evident during early season events when the 
Sea-surface temperature decreases while the net surface heat flux is down- 
ward. This means that a forecast scheme must be capable of predicting 

The changes in the mixed layer depth to predict the sea-surface tempera- 
Ture evolution. 

4, The data sets examined in this study suggest that at OWS P the 
large forcing events are largely dominated by mechanical mixing. This 
is evidenced by the entrainment heat ae exceeding the surface heat 
flux, and the increase of potential energy of the upper ocean. At OWS V 
and OWS N, however, the large majority of the cases showed that the 
strong forcing events are dominated by the convective process. 

3. The adequacy with which the EFT model simulated the evolution 
during these strong events suggests that a properly parameterized turbu- 
lent bulk model will be useful for predicting the thermal structure 
changes over a period of a few weeks. Even during periods when non-local! 
processes were important, the EFT model was capable of predicting changes 
in the thermal gradient. This property should be useful for estimating 
the changing acoustic properties of the upper ocean due to storm activity. 

6. In agreement with the findings of Dorman (1974) the modeling of 
data sets at OWS N and OWS V suggests that the inclusion of salinity 
effects may not be necessary to simulate the upper ocean thermal struc- 
ture changes in the subtropics. Moreover, the results at OWS P suggest 


that a model that neglects salinity effects is capable of simulating the 
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changes in the temperature structure during the early cooling season 
in the subarctic region as well. 

7. The examination of the seasonal trends in the mixed layer evolu- 
tion suggests that a large percentage of the changes are understandable 
in terms of one-dimensional processes. Additionally, the capability of 
the modified EFT model to simulate layer retreat suggests that it should 
be useful for simulating thermal structure changes during the spring and 
Summer heating seasons. 

8. The examination of the forcing terms showed that the distribu- 
tion of the marine winds and turbulent heat fluxes are non-Gaussian at 
the three ocean weather stations. The analysis showed that although the 
magnitude of the mechanical forcing is largest at OWS P and smallest at 
OWS N, its frequency distribution is similar at all stations. The dis- 
tribution of the turbulent heat fluxes were similar at OWS P and OWS V, 
but different at OWS N. It was demonstrated that this is related to the 
fact that the large forcing events at OWS P and OWS V are very similar 
and are closely correlated to the extratropical cyclones that pass in 
close proximity to these two stations. At OWS N, however, cyclone activ- 
ity is rare and the events are related to a pulsing of the mean flow. 
This pulsing is probably due to a strengthening of the north-south pres- 
sure gradient as the large storms pass north of QWS N. 

faa ine Tinal conciusicen 1s that The Timing ef the large fall vamd 
winter events and the characteristics of the underlying thermal structure 
are important to the seasonal evolution of the upper ocean. It was 
demonstrated that the strong forcing events that occur early in the 
cooling season result in a much larger mixed layer response than compar- 


able events occurring later in the season. This is because less turbulent 
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kinetic energy is available for entrainment when the layer is deep (dis- 
Sipation enhancement). It was further demonstrated that changes in the 
mixed layer will be quite different depending on the strength of the 
thermocline. Cursory examination of records from OWS P suggests that 
the strength or weakness of the seasonal thermocline depends on the 
characteristics of the forcing during the spring and summer heating sea- 
son. It is possible, therefore, that the large forcing events are 
equally important in the formation of the thermocline, during the spring 
and summer, as they are in its subsequent erosion during the fall and 
winter. It is recommended that future research be conducted during the 
spring and summer periods to demonstrate the role of strong events during 


These seasons. 
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APPENDIX A 


COMPUTATIONAL FORMULAS FOR SURFACE FORCING 


Seen DIATIVE FLUXES AT THE SEA SURFACE 

The heat gained or lost by the oceans at the air-sea interface by 
radiant energy falls into two spectral regions. The first, ranging 
from O.| to 4 microns oe cm), is the short wave radiation received 
from the sun. The long wave, 4 to 50 microns, is commonly known as the 
back radiation and represents a net loss of heat from the ocean to the 
atmosphere or space. These two spectral ranges are virtually exclusive, 
thus permitting computations of these radiative fluxes to be performed 
separately. 

A number of empirical formulas are available for computing the in- 
solation arriving at the sea surface and Reed (1975) has reviewed and 
evaluated the most commonly used expressions. Reed found that results 
obtained from the formula developed by Seckel and Beaudry (1973), were 
consistently in better agreement with data collected at five coastal] 
stations. Therefore this expression is used to calculate the clear-sky 


radiation Oy ine Ti iSesauidy. 


Q5 = nage A SOSie+ 8 Sind + A, COS Zoe 5. Sin 20 (A-1) 
where oF is in langley (ly) per day, 
b = 2. Galen 


and Tt is the julian day of the year. The coefficients CAS» etc.) were 


calculated by a harmonic representation of the values presented in the 
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Smithsonian Meteorological Tables (List, 1958). This clear-sky value 
must be corrected for the presence of clouds and reflection from the 
sea surface. Gunter Seckel (personal communication) has suggested that 
the cubic cloud correction of Laevastu (1/960) and a reflection coeffi- 
cient modeled after Anderson (1952) are suitable at the ocean stations. 


Therefore, 


- = 6 K Cl = R) (A-2) 


calculates the solar energy penetrating the air-sea interface and 


(=i aenoomee 


The coefficient C is the total observed cloud cover (in tenths) and 
a is the mid-day elevation angle of the sun. The constants, a and 
Eeeeecre adopted from fabata (1964), and for C< 0.5 , a = 0.35 and 
Meeeno-42 , while for C > 0.5 , a on Ore 4 |AMaes) ((onmemoneme dO sc) 

The largest source of error in (A=-2) is the parameterization of the 
effects of cloud cover and the subjectivity involved in observation. 
Moreover, the application of this expression for averages less than mean 
monthly values introduces another possible source of error, as discussed 
by Reed and Halpern (1975). In this research (A-2) was used to estimate 
o ona daily basis with C taken as the mean cloud cover during the 
daytime. 

The net long wave radiation (Q,? is a function of the radiation 
emitted from the sea surface to the atmosphere, minus the energy radiated 
from the air mass and absorbed by the ocean. Both of these quantities 


depend upon the fourth power of the absolute temperature of the emitting 





body (Stephan Boltzman Law) with suitable correction factors for cloud 
cover and vapor content of the atmosphere. A representative formula 


reported by Husby and Seckel (1975) is 


a 2 


Q, = 1.14x10° (273.164) (0.39-0.05 vE_)(1-0.6C*) (A-3) 
where 4 is in ly/day, Ve is the sea-surface temperature (°C), and 
Bi is the saturated vapor pressure of the atmosphere (at a height of 

lO m) in millibars. This vapor pressure was calculated using the Goff- 
Gratch (1946) formulation of the ClausiusClapeyron equation, using the 
dew point temperature (Ty) as the entering argument. Equation (A=-3) is 
the modified Brunt (1932) formula with the empirical constants of Budyko 


(1956), and the largest uncertainties are introduced through the cloud 


correction factor and the use of overland constants. 


B. TURBULENT FLUXES OF HEAT AND MOMENTUM 
The turbulent fluxes of latent heat (Q,)> sensible heat (Q,), and 
momentum ee at the air-sea interface were represented by the so-called 


bulk aerodynamic formulas. 


& Z Z 
ie = omen Cu. x 10°) (dynes/cm™ ) (A-4) 
on = 3,767 Cy (O57. EE.) u. (ly/day) (A-5) 
Q, = 2,488 C, (T.-T,) us (ly/day) G26) 


where u, is the mean wind speed (m/sec), i is the air temperature 
ec), E is the saturated vapor pressure of the marine air directly in 
contact with the sea surface (0.98 corrects for salt effects), and Cy 


is the non-dimensional! drag coefficient. A constant drag coefficient 
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igi. X on is used in all computations in this thesis and is consis- 
tent with the range of values reported in the literature. 

The accuracy of these expressions has been the subject of many de- 
tailed studies and the main sources of error are the underlying assump- 
tions of a neutrally stable atmosphere and constant and equal exchange 
coefficients (moisture, momentum, and heat). Businger, et al. (1971), 
macm overland values, and Paulson, et al. (19/2), from data collected 
at sea, have demonstrated that the moisture and heat coefficients are 
nearly equal but quite different from the coefficients of momentum ex- 
change. Furthermore, studies by Deardorff (1968), DeLeonibus (1971), 
and Davidson (1974) have found that these coefficients are very depen- 
dent on the stability of the marine boundary layer and the roughness of 
The sea surface. Nevertheless, they afford the only practical means 
for computing these fluxes, using the meteorological observations avail- 
able in the ocean weather station file, and are used throughout this 


research. 


io 





APPENDIX B 
THE NUMERICAL SCHEME FOR THE ONE- 
DIMENSIONAL TURBULENT BULK MODELS 

A simple numerical scheme was developed which is capable of incor- 
porating the various assumptions of the Kraus-Turner (1967), Elsberry, 
ef al. (1976), and Kim (1976) models in a consistent manner. This rou- 
tine is a modification of the algorithm presented by Thompson (1976), 
and is designed to calculate the potential energy modifications due to 
the vertical fluxes of heat and turbulent kinetic energy. 

The NODC mechanical bathythermographs that were used to initialize 
and validate the models were digitized in five meter increments starting 
at the sea~surface. In the model, this temperature profile is stored in 
N equally-spaced grid intervals [-(n-1) AZ, -nAZ], where n=1,2,...,N, 
and AZ = 2.5m. Fora profile in which The temperature is wel! mtxed 
to a depth of -nAZ , the first n-values of JT are equal to qT eee tee 
Ing initialization it was assumed that The temperature varied linearly 
with depth in each five meter increment in the BT profile, and the mode! 
profile was chosen to conserve the heat content of the BT profile. 


In the model the potential energy (PE) per unit area may be expressed 


as 


O 
Ee — “P90 f HZdz (B-1) 
-D 


where D 1s a depth which is deeper than the maximum penetration of the 
vertical turbulent processes (typically the deepest level in the model). 
If non-local processes are smal! relative to the vertical fluxes, and 


if density changes due to salinity changes can be ignored, and a is 
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constant, changes in the potential energy calculated by (B-!) will be 
representative of the changes in the potential energy per unit area 
of the ocean. 

After mixing the top n grid intervals a further mixing of the layer 
to a depth -(ntl)AZ will result in a change in potential energy (APE) 


of the column by 


ier! Z 
APE(n}) = x PgantAZ) oP Tate (B-2) 
Zor Tes - ae > 0, this mixing increases the potential energy and 
APE (n) will represent the amount of turbulent kinetic energy that must 


be expended to accomplish the mixing. However, if the column is unstable 
ol. - eer < 0), this mixing will release potential energy, and APE (n) 
will represent the turbulent kinetic energy generated by free convection. 


At the beginning of each time step (At = | hr) The surface heat 


fluxes are added to obtain 


eT +AT ) 


* * < 
! T (440400, (4 40560, 7 1-9, (AZ,t 1/06 (B-3) 


p 


* % 
T (4) +0419. (-nAZ, t )-Q.(- (nt! AZ, 4 ) 1/9 ,C (Baa) 


eT taTt) 
n p 


ie@em=2Z,2,~«..,N, and + =f + = At . The effective insolation (Q.? was 
distributed by assuming 50% absorption in the first meter and The re- 

mainder taken to decay as exp(-yZ). The average extinction coefficient 
(y) was assumed constant at 0.005 ae which is the value used by Den- 
man and Miyake (1973). If the resulting temperavure ororile is unstable 


(T < T5) the first two intervals are mixed and the resulting turbulent 


kinetic energy (APE) generated by the convection is calculated according 


ie 5 


“be 





to (B-!). This process continues until a stable temperature profile is 


established (Ty el ) and the upper layer is isothermal to a depth 


k+ | 
Z= -kAZ . At the end of this evolution the potential energy of the 


column will be reduced by 
k- | 
APE _(n) ; (B=) 
n= | 


which is also equal to the total generation of turbulent kinetic energy 
by free convection. 

Further mixing will require an expenditure of turbulent kinetic 
energy. 

The algorithm must also be modified to account for dissipation en- 
hancement (EFT and KIM), the storage of turbulent kinetic energy (KIM), 
and the fraction of convectively generated turbulent kinetic energy that 
is utilized for entrainment (KT, EFT, and KIM). We start by defining 


the total amount of turbulent Kinetic energy, E.(n), available for mix- 


- 
ing the first n levels with level n+l as 
Ee Gi) S= 7b ue eater nD (= iene (B-6) 
T m eS p S 
The mechanically generated turbulent kinetic energy is EM) : ES is 


The fraction of (B=-5) available for entrainment, Sey is the amount 
of turbulent kinetic energy previously expended to mix the layer to 
level n, and E(n) is the amount of turbulent kinetic energy that is 
stored as the layer deepens to a depth -nAZ . 


For the KIM model, 


a a 
emt) = (Cis: OW x (tT) - P,D,,nAZ) ie (B-7) 
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where D., is a constant background dissipation equal to 2 x 1074 aie 


poet. For the KIT and EFT models, 


= 
Et) Co Wy (+t ) exp(-nAZ/Z] At , (B=-8) 


with 


co for KT 
Ze (B-9) 

DOsmee TOmeer |: 
Therefore (B-7), (B-8), and (B=9) express the depth dependence of dissi- 
pation (dissipation enhancement) formulated by Kim (1976) and Elsberry, 
eee. (1976). In the Kraus and Turner (1967) model dissipation is 

* 

neglected, and the amount of surface production, DWy Ct ) , available 


for mixing is independent of depth. 


Further, 
k= | 
—E = 2s ABE (Cn? (B=-10) 
e aa e 


where r is the fraction of the convectively generated turbulent kinetic 


energy utilized for entrainment. In all models used in this study fr = 


-15 following Gill and Turner (1976). Therefore the KT model here is 
equivalent to a modified version discussed by Gill and Turner (1976). 
Finally, 
n= | 
an) * Deane (B-11) 
i=k 


67 


vy 


i 





and | 


n= 
3. * 
> 4.5 0 AZ Wome emote e114 
i=2 
a= cea) 


O , for KT and EFT 


Thus (B-I1) is used to calculate the turbulent kinetic energy expended 
to deepen the layer from the free-convection depth to -nAZ . For KIM, 
(B-12) calculates the amount of turbulent kinetic energy stored as the 
mixed layer deepens beyond the mixed layer depth of the previous Time 
Paeor i.e., SAZ(t+AT) > h(t+) J]. 


nee ENS APE (n) , there is enough energy to mix T through ag 


T | ° 


lf there is insufficient energy to mix completely, E-(n) < APE (n) . 


then we follow Thompson (1976) and partially mix; I[.e., set 


Messe 
ee 
(peel ve Sm) 
m 
T= Lat, + (n-a) TJ/n 
Tl= me 
T a ines Cl —ay Ta 
Siepeset 1. = 1, for i=l,2,...n, 1 = T' . The mixed-layer tempera- 


| ice! 


Ture at each time step is equal to T, but since the algorithm does not 
calculate the mixed layer depth, it is defined as the depth at which the 
temperature is 0.2°C less than T | - 
The relative Importance of mechanical mixing and free convection to 
the evolution of the mixed layer were investigated by comparing the rela- 
mive contribution of APE and APE. to APE. An important part of 


This study was to quantify the relative magnitudes of the surface and 


entrainment heat fluxes. The surface heat fluxes were calculated, as 


168 


pi 


\e 





described in Appendix A, from the surface marine observations. fhe en- 
trainment heat flux (A ae AT) was calculated, in the models, as the 
heat gained by the profile (during each time step) in fevels n+l to ntm 


as the mixed layer deepens between levels n and ntm. 
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